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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

\\\ 

The  research  described  in  this  report  was  undertaken  to  extend  our 
knowledge  of  the  chemical  transformations  of  hydrazine  vapor  in  the 
environment.  It  is  known  that  both  gas-phase  and  surface-catalyzed  reactions 
can  occur  with  this  fuel.^In  the  immediate  vicinity  of  the  ground,  other 
surface  processes,  such  as  adsorption,  can  also  be  important. 

^To  determine  the  eventual  fate  of  hydrazine  vapors  which  escape  into  the 
atmosphere  through  normal  operations  or  during  accidental  liquid  spills,  all 
of  these  processes  must  be  investigated,  This  study  was  undertaken  in  a 
chamber  designed  to  simulate  atmospheric  conditions,  but  allowing  for  some 
control  over  important  variables  to  address  some  of  the  potential  chemical 
transformations. 


B.  BACKGROUND 


\iydrazine  is  used  extensively  in  small  thrusters  for  satellites,  as  a 
fuel  cell  reactant,  in  emergency  power  generating  units,  and  as  a  liquid 
rocket  fuel  component.  Both  routine  handling  operations  and  accidental 
spills  can  produce  fuel  vapors  which  constitute  a  substantial  health  risk* 
Hydrazine  is  classified  as  an  animal  carcinogen  (Reference  1)  and  a  suspect 
human  carcinogen .  <Referenc« -2)^  This  toxicity  is  reflected  in  its  low 
threshold  limit  value,  0.1  parts  per  million  (ppm)  (Reference  3),  and  has 
prompted  numerous  studies  of  the  fate  of  hydrazine  in  both  atmospheric  and 
terrestrial  environments.  /jf,*-  t  fC  f  WT-r.y,  73?>rtC 


The  autoxldation  of  hydrazine  vapor  has  been  the  subject  of  several 
previous  studies.  Bowen  and  Birley  investigated  what  they  believed  to  be  the 
vapor-phase  reaction  between  hydrazine  and  oxygen  in  a  Pyrex  reaction  vessel 
(Reference  4).  They  concluded  that  from  100  °C  to  160  °C  the  reaction 
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was  first  order  in  both  hydrazine  and  oxygen  and  that  it  proceeded  by  a  five- 
or  six-step  successive  chain  reaction.  The  overall  reaction  stoichiometry 
was  shown  to  be: 


N2H4  +  °2  - ►  N2  +  2H2°  (1) 

They  further  noted  that  the  rate  of  the  reaction  was  a  function  of  the 
surface  area  of  the  reaction  vessel. 

Winning  (Reference  5)  conducted  a  similar  study  in  2-L  and  0.5-L  quarts 
reaction  vessels  from  20  °C  to  185  °C.  He  also  found  that  the  rate  of 
the  reaction  was  directly  proportional  to  the  surface  area  of  the  reaction 
vessel.  His  results  also  showed  that  the  reaction  was  first- order  with 

A 

respect  to  both  hydrazine  and  oxygen.  He  determined  that  neither  added 
nitrogen  or  water  vapor  had  any  effect  on  the  reaction  rate.  However,  he 
interpreted  his  results  by  invoking  a  nonchain  bimolecuiar  surface  reaction. 
He  also  noted  the  production  of  some  ammonia  at  temperatures  above  115  °C. 

Stone  conducted  two  studies  of  the  aut oxidation  of  hydrazine  vapor  in 
variou8“sized  Pyrex  containers  (References  6  and  7).  All  experiments  were 
conducted  at  25  °C.  These  studies  gave  a  reaction  stoichiometry  in 
agreement  with  equation  (1).  They  also  showed  the  dependence  of  reaction 
rate  on  aurface-to-volume  ratio  and  on  the  surface  composition.  Small 
amounts  of  ammonia  were  produced  in  most  experimental  runs. 

Bellerby  (References  8  and  9)  analyzed  the  studies  of  earlier  researchers 
and  deduced  a  mechanism  which  accounts  for  their  major  experimental  observa¬ 
tions.  He  proposed  that  the  initial  formation  of  a  hydroperoxide  by  a  chain 
process  was  followed  by  wall  reaction  to  give  the  observed  products. 

An  entire  doctoral  study,  "The  Vapor  Phase  Oxidation  of  Hydrazine, M  was 
recently  completed  by  Moody  (Reference  10).  He  studied  the  effects  of 
temperature,  surface-to-volume  ratio,  total  reaction  pressure,  and  the 
composition  of  the  surface  of  the  reaction  vessel.  The  reaction  vessels  used 
in  this  work  were  relatively  small,  0.3  -  1.9  L,  and  constructed  from  Pyrex. 


In  these  vessels,  Moody  found  that  there  was  an  unusual  negative  temperature 
coefficient  in  the  reaction  rate.  He  concluded  that  the  surfaces  of  the 
reaction  vessels  play  a  dominant  role  and  that  there  is  ample  evidence  for  a 
free-radical  chain  reaction  mechanism. 

Pitts  and  co-workers  conducted  two  studies  of  the  reactions  of  hydrazines 
under  simulated  urban  tropospheric  conditions  (References  11  and  12).  These 
studies  were  conducted  in  large  Teflon  film  reaction  chambers  (volumes  were 
30,000,  6,400,  and  3,800  liters)  at  room  temperature.  The  reactions  of 
hydrazine  and  various  derivatives  with  hydroxyl  radicals,  oxides  of  nitrogen, 
ozone,  nitric  acid,  and  formaldehyde  were  studied.  It  was  determined  that 
hydrazine  has  an  expected  gas-phase  lifetime  of  about  30  minutes  under 
typical  urban  atmospheric  conditions.  These  investigators  also  used  radical 
traps  and  chemical  tracers  to  determine  details  of  the  reaction  mechanisms 
involved  in  the  atmospheric  oxidation  of  hydrazines. 

Naik  and  co-workers  also  conducted  a  study  of  hydrazine  air  oxidation  in 
a  Teflon  film  reaction  chamber  (References  13  and  14).  The  chamber  they  used 
was  smaller  (about  320  L)  and  their  studies  centered  on  the  effects  of  added 
surfaces  on  the  rate  of  hydrazine  oxidation  in  air.  They  determined  that 
among  the  added  metal  surfaces  studied,  copper  was  by  far  the  most  reactive, 
followed  by  painted  aluminum,  aluminum  foil,  stainless  steel,  and  aluminum 
plate.  They  also  noted  that  the  half-life  of  the  hydrazine  was  a  strong 
function  of  the  relative  humidity  of  the  air  in  the  chamber. 

Two  additional  studies  have  been  conducted  in  smaller  (-100-L)  synthetic 
film  bags.  The  first  was  conducted  by  Lin  (Reference  15),  who  determined  the 
decay  rates  of  hydrazine  in  nitrogen  and  in  dry  air  in  Teflon  and  Tedlar 
bags.  The  half-life  of  the  hydrazine  was  notably  shorter  in  the  Tedlar  bag. 
In  a  follow-up  to  some  of  his  earlier  work,  Naik  studied  the  decay  of 
hydrazine  in  a  Teflon  bag  (Reference  16).  In  this  study,  he  determined  the 
half-life  of  hydrazine  in  dry  air  as  a  function  of  the  surface-to-volume 
ratio  of  the  Teflon  bag.  By  extrapolating  to  a  surface-to-volume  ratio  of 
zero,  he  estimated  that  the  homogeneous  half-life  of  hydrazine  in  air  is 
about  60  hours. 
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Previous  studies  have  shown  that  the  relative  magnitudes  of  Msnidity, 
surface  adsorption,  permeation  through  walls,  surface-to-volume  ratio  and 
actual  oxidation  as  factors  affecting  the  decay  of  hydrazine  In  Teflon  film 
reaction  chambers  are  still  questionable.  This  study  was  undertaken  to 
further  investigate  these  effects  with  the  aim  of  attempting  to  clarify  the 
relative  importance  of  these  various  loss  processes. 

C.  SCOPE 

The  study  included  a  thorough  examination  of  the  behavior  of  hydrazine 
vapor  in  a  Teflon  film  reaction  chamber.  It  included  the  design  and 
construction  of  the  chamber  and  a  determination  of  its  optimum  operating 
conditions.  It  also  included  the  development  of  a  comprehensive,  kinetic 
model  to  account  for  all  potential  molecular  transformations,  both  chemical 
and  physical.  Finally,  it  included  chamber  studies  of  the  decay  of  hydrazine 
vapor  under  four  types  of  experimental  conditions:  dry  and  humid  nitrogen, 
and  dry  and  humid  air. 


SECTION  II 


EXPERIMENTAL  PROCEDURES 


A.  REACTION  CHAMBER 

The  reaction  chamber  was  constructed  by  heat-sealing  pieces  of  DuPont  FEP 
Teflon  film.  The  film  thickness  was  .005  inches  and  the  resulting  chamber 
dimensions  were  approximately  0.87  by  0.66  by  0.56  meters,  giving  a  volume  of 
about  320  liters.  The  chamber  was  assembled  by  technicians  at  the  NASA  White 
Sands  Test  Facility,  White  Sands,  NM.  One  of  the  longer  sides  of  the  bag  had 
a  0.3-  by  0.3-meter  access  port  which  was  constructed  by  the  same 
heat-sealing  technique.  This  port  extended  out  from  the  side  of  the  chamber 
about  0.5  meter  in  the  form  of  a  square  sleeve.  The  access  port  was  sealed 
during  chamber  experiments  by  carefully  folding  the  Teflon  sheets  and  tightly 
clamping  them  between  two  wooden  supports.  The  comers  of  the  chamber  were 
reinforced  with  strips  of  2-inch  Mylar  tape  (3M  core  series  2-0300,  3M  Co., 
St.  Paul,  MN),  including  the  access  port.  The  chamber  is  shown  schematically 
in  Figure  1. 

The  chamber  was  supported  by  securing  its  two  longer  upper  edges  to 
0.5-inch  aluminum  rods  with  Mylar  tape.  The  aluminum  rods  were  then 
suspended  with  wire  from  structural  members  of  a  supporting  rectangular  frame 
which  measured  0.75  meters  by  1.18  meters  by  1.37  meters.  The  frame  members 
were  hollow,  1.0- inch  square  aluminum  bar  stock  welded  together.  The 
aluminum  frame  was  covered  with  1/8-inch  Plexiglas  and  placed  on  a  wooden 
table,  completely  enclosing  the  Teflon  chamber.  A  1/6  hp  motor  (Lightning 
Model  L,  1550  B/min)  was  mounted  on  the  structural  aluminum  frame  members 
above  the  Teflon  chamber.  The  shaft  of  the  motor  was  connected  to  a  6-inch, 
Teflon-coated  fan  blade  (located  inside  the  Teflon  chamber)  with  a  sealed, 
rotatable  feed  through.  The  Plexiglas  box  was  fitted  with  four  1/4-inch 
Nylon  bulkhead  fittings  with  1/4-inch  PVC  tubing  leading  to  the  system  gas 
distribution  manifold.  These  ports  were  used  to  direct  gas  flow  into  the 
annular  space  between  the  Plexiglas  box  and  the  Teflon  chamber.  One  end  of 
the  Plexiglas  was  hinged  for  access  to  the  optical  components.  The  other  end 
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Figure  1.  Schematic  Drawing  of  the  Teflon  Chamber . 


of  the  box  had  a  15  by  40  cm  hinged  opening  for  both  beam  path  and  path 
length  adjustments.  There  was  an  additional  hinged  opening,  40  by  40  cm  for 
access  to  the  Teflon  chamber  sample  port.  There  were  other  openings  for  the 
infrared  beam,  electrical  cords,  gas  lines,  and  the  exhaust  port.  The 
Plexiglas-covered  aluminum  frame  is  shown  in  Figure  2. 

The  chamber  had  three  other  access  ports.  One  was  a  standard  24/40 
ground  glass  joint  located  in  the  lower  comer  of  the  chamber  and  extended 
through  the  Plexiglas  shell.  This  port  was  used  for  inserting  a  thermistor 
probe  for  temperature  measurements  of  the  chamber  contents.  The  temperature 
inside  the  chamber  was  continuously  monitored  by  recording  the  signal  from  a 
thermistor  (Yellow  Springs  Instruments  Model  401)  thermometer  (Yellow  Springs 
Instruments  Model  44TD)  with  a  strip  chart  recorder  (Varian  Model  9176).  The 
recorder  response  was  calibrated  using  known  water  bath  temperatures. 

The  second  port  was  a  0.5-inch  glass  tube,  located  in  the  upper  surface 
of  the  chamber  and  connected  to  a  butterfly  valve  on  the  top  Plexiglas 
surface  of  the  chamber  cover.  This  valve  was  connected  to  the  PVC  fume  hood 
extension  through  a  2-inch  piece  of  Tygon  tubing.  The  third  port,  the  sample 
inlet  port,  was  a  piece  of  10  mm  0D  Pyrex  tubing  which  entered  the  Teflon 
chamber  from  the  side  opposite  the  temperature,  probe  port.  Each  of  these 
glass  ports  was  Inserted  through  the  chamber  wall  and  sealed  with  carefully 
placed  strips  of  Mylar  tape. 

The  chamber  was  checked  for  leaks  by  first  being  flushed  with  helium  gas, 
then  inflated  to  slightly  over  its  equilibrium  "reasure.  A  gas  leak  detector 
(Gow  Mac  Model  21-150)  was  then  used  to  inspect  the  seams  and  corners  for 
leaking  helium.  A  number  of  leaks  found  in  the  chamber  were  repaired  by 
pressing  small  strips  of  Mylar  tape  over  the  leaking  areas.  At  some  of  the 
comer  Joints,  silicon  glue  was  squeezed  into  an  area  whose  geometry  did  not 
allow  the  efficient  use  of  tape.  Later,  pressure-sensitive  tape  made  from 
unsupported  Teflon  TFE  film  (Berghof /America,  Inc.,  Type  15078)  was  used  to 
repair  leaks.  This  tape  was  more  flexible  and  provided  better  seels  than  the 
M/lar  tape. 
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Figure  2.  Schematic  Drawing  of  the  Plexiglas-Covered  Aluminum  Support  Structure  for  the 
Teflon  Chamber. 


Also  contained  within  the  Plexiglas-covered  aluminum  frame  was  the  sample 
inlet  manifold.  This  manifold  consisted  of  three  Kontes  10  mm  high  vacuum 
valves  with  Teflon  seals  for  admitting  gaseous  samples  and  a  silicon  septum 
inlet  for  injecting  liquid  or  gas  samples.  The  manifold  was  connected  to  a 
500  mL  mixing  bulb  and  then  proceeded  through  10  mm  Pyrex  tubing  to  the 
Teflon  chamber.  The  sample  manifold  was  heated  with  heating  tape  controlled 
by  variable  autotransformers.  The  sample  manifold  is  shown  in  Figure  3. 

B.  PURE  AIR  SYSTEM 

Pure,  dry  air  for  chamber  purging  and  filling  was  generated  by  a  locally 
constructed  system.  A  7.5  hp  air  compressor  (Ingersoll-Rand  Model  T30)  was 
used  to  deliver  compressed  air  to  a  noncycling  refrigerated  air  dryer 
(Zeks-Therm  Model  35NCA) .  This  unit  delivered  air  at  a  dew  point  of  35  UF 
to  a  catalytic  air  purifier  (Aadco  Model  737-15A).  The  air-purifier  unit  was 
capable  of  delivering  250  L/min  of  air  with  a  less  than  one  part-per-billion 
ozone,  methane,  hydrocarbons,  carbon  monoxide,  hydrogen  sulfide,  sulfur 
dioxide,  or  fluorocarbons  at  a  dewpoint  of  -60  °F.  The  pure  air  was 
delivered  to  the  chamber  through  1/2-inch  PVC  tubing  (Imperial  Eastman 
"Poly-Flo'4'  tubing,  type  88-P-1/2).  The  pure  air  system  and  distribution 
manifold  are  shown  schematically  in  Figure  4. 

C.  FT-IR  SPECTROMETER  AND  LONG-PATH  CHAMBER  OPTICS 

Molecular  species  of  interest  in  this  study  were  monitored  by  Fourier 
transform  infrared  (FT-IR)  spectroscopy  using  a  Piaentel-type  (Reference  17) 
multiple-reflection  mirror  system.  The  principal  mirrors  were  located  within 
the  Teflon  film  chamber.  These  mirrors  consisted  of  an  in-focus  (nesting) 
mirror  cut  from  an  8-inch  spherical  mirror  and  two  4-inch  out-of-focus 
(collecting)  spherical  mirrors,  all  with  a  common  radius  of  curvature  of  31.5 
inches.  These  mirrors  were  standard  1/4-wave,  first-surface  Pyrex  mirrors 
aluminized  and  overcoated  with  silicon  monoxide  (A.  Jaegers  Co.).  They  were 
used  as  received.  The  corner  reflector  was  made  by  cutting  a  section  of  the 
nesting  mirror  and  gluing  it  at  right  angles  to  a  section  of  flat  mirror  with 
epoxy.  The  corner  reflector  was  glued  on  a  standard  optical  post  assembly 
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Schematic  Drawing  of  the  Chamber  Sample  Injection  and  Introduction  System. 
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Figure  4.  Schematic  Drawing  of  the  Laboratory  Pure  Air  Generation  and  Delivery  System. 


for  rotational  and  height  adjustability.  The  other  mirrors  were  glued  to 
positioning  mounts  (Newport  Model  MM-2)  fastened  to  standard  optical  posts, 
enabling  the  required  x-y-z  adjustments. 

The  total  infrared  beam  path  within  the  chamber  was  determined  by 
counting  the  red  image  dots  from  a  helium-neon  laser  whose  optical  path  was 
coincident  with  that  of  the  infrared  beam  from  the  FT-IR  system.  The 
long-path  optical  system  was  operated  at  46  passes,  giving  a  path  length  of 
36.8  meters. 

The  incoming  infrared  beam  was  directed  into  the  chamber  and  then  onto  a 
mercury-cadmium-telluride  detector,  operated  at  77  K  by  a  series  of  external 
mirrors,  as  shown  in  Figure  5.  The  FT-IR  spectrometer  employed  was  a  Nicolet 
Model  6000C.  This  spectrometer  has  a  maximum  resolution  of  0.06  cm  \  and 
was  upgraded  to  include  a  raster  scan  display,  24-Mbyte  hard  disk,  1-Mbyte, 
8-inch  floppy  disk,  digital  plotter,  and  RS-232  dot  matrix  printer.  Thus*, 
configured,  it  operated  on  the  more  recent  NXCOS  software. 

D.  MATERIALS 

The  hydrazine  was  anhydrous,  fuel  grade  (98  percent  pure)  from  Rocky 
Mountain  Arsenal.  It  was  used  as  received.  Methane  (CH^)  gas  was  Fisher 
high-purity  grade  (99.99  percent)  and  sulfur  hexafluoride  (SF^)  was  Air 
Products  instrument  grade  (99.99  percent). 

E.  METHODS  AND  PROCEDURE 

1.  Sample  Introduction 

To  conduct  an  experiment  in  dry  air,  the  following  procedure  was 
used.  The  annular  apace  between  the  Plexiglas  box  and  the  Teflon  film  bag 
was  purged  with  dry,  purified  air  from  the  pure  air  system  at  20  L/min  for  at 
least  60  minutes.  At  the  same  time,  the  Teflon  film  chamber  was  purged  with 
a  separate,  10  L/min  flow  of  pure,  dry  air.  The  chamber  mixing  fan  was 
operated  during  the  entive  purging  process.  The  flow  of  air  into  the  chamber 
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was  manually  balanced  with  the  removal  of  aiF  from  the  chamber  via  the 
exhaust  port.  After  60  to  90  minutes  of  purging,  the  humidity  in  the  annular 
space  of  the  box  and  the  humidity  in  the  chamber  were  measured  electronically 
with  a  thin  film  capacitance  device  (Vaisala  Model  HMI  31  humidity  and  Model 
HMP  31  UT  probe) .  The  readings  were  always  less  than  2  percent  relative 
humidity.  The  chamber  was  now  ready  for  sample  introduction. 

To  inject  a  sample,  the  purge  flow  in  the  annular  space  was  lowered 
to  about  10  L/min  to  maintain  the  low  humidity,  and  flow  to  the  chamber  was 
stopped.  The  volume  of  the  chamber  was  lowered  by  about  20  percent  by 
opening  the  exhaust  port.  At  this  point,  a  reference  spectrum  was  taken  with 
the  FT-IR  system  (see  the  following  section  for  details  on  how  the  FT-IR 
spectrometer  was  operated).  Both  gas  and  liquid  samples  were  injected  with 
syringes  through  a  silicon  septum  in  a  cylindrical,  flow-through  Pyrex 
sampling  vessel.  The  sampling  vessel  was  heated  to  about  50  C  with  heat 
tape  for  hydrazine  injections.  Samples  were  flushed  into  the  chamber  with 
pure  air  flowing  at  about  2  L/min.  The  mixing  fan  was  run  continuously 
during  sample  injections.  Pure  air  flow  was  maintained  until  the  chamber  was 
reinflated.  Care  was  taken  to  see  that  the  chamber  was  actually  somewhat 
underinflated  to  prevent  any  physical  force  being  exerted  by  the  walls  of  the 
chamber  which  might  have  unnecessarily  promoted  various  loss  processes. 

To  conduct  an  experiment  J.n  humid  air,  the  procedure  outlined  above 
was  followed,  except  that  the  pure  air  flows  for  both  the  annular  box  space 
and  the  chamber  were  routed  through  a  10  cm  diameter  by  75  cm  high  Plexiglas 
column  filled  with  Norton  16  mm,  polypropylene  Pall  rings,  and  distilled 
water.  Flow  through  this  column  resulted  in  relative  humidity  levels  of 
about  60  percent  in  the  chamber  and  box  system. 

Experiments  which  used  dry  nitrogen  also  used  the  same  procedure 
outlined  in  the  beginning  of  this  section.  The  nitrogen  was  obtained  from 
the  boiloff  of  liquid  nitrogen. 


2.  FT-IR  Spectrometer  Operation 


Molecular  species  of  interest  in  this  study  were  monitored  by  record¬ 
ing  their  absorption  spectra  as  a  function  of  time  with  FT-IR  spectroscopy. 
A  pathlength  of  36.8  meters  was  set  in  the  Teflon  chamber  and  a  resolution  of 
1.0  cm  ^  was  employed.  Data  were  recorded  in  the  mid- infrared  from  4000  - 
700  cm  *  using  a  mercury-cadmium-telluride  (MCT)  photoconductive  detector 
(Infrared  Associates  Model  G-2229). 

A  background  spectrum  for  a  particular  experiment  was  obtained  by 
co-adding  256  instrument  scans  and  then  storing  the  result.  Immediately 
after  the  injection  of  a  sample,  the  instrument  was  placed  under  software 
control  for  automatic  data  acquisition,  processing,  storage,  and  plotting. 
This  gave  a  spectrum  within  2  to  5  minutes  of  sample  injection.  Three 
additional  spectra  were  acquired  in  the  same  manner  at  10-minute  intervals. 
The  instrument  was  then  directed  to  increase  the  data  acquisition  interval  to 
a  longer,  operator-controlled  value,  and  continue  to  acquire,  process,  store, 
and  plot  data  for  up  to  100  additional  intervals.  Generally,  these  longer 
intervals  were  30  to  60  minutes,  using  60  instrument  scans.  Experiments  were 
conducted  for  periods  of  12  to  60  hours. 

The  spectral  acquisition  software,  written  using  a  series  of  Nicolet 
MACRO  commands,  produced  single-beam  spectra  stored  in  successive  files  in 
the  “scratch"  area  of  the  system  hard  disk.  These  single-beam  spectra  were 
ratloed  against  the  background  spectrum,  converted  to  absorbance  mode,  and 
analyzed  to  determine  the  absorption  at  specified  analytical  wavelengths  for 
each  molecular  species  of  Interest.  Finally,  the  absorption  values  were 
printed  out,  along  with  time  and  wavelength  information.  All  these 
operations  were  handled  automatically  by  the  instrument's  MACRO  software. 
(See  Appendix  A  for  program  listings). 

The  automation  of  data  acquisition  enabled  long  experiments  to  be 
conducted  without  constant  operator  attention.  However,  this  automation  also 
caused  a  false  temperature  dependence  to  be  observed  In  the  data.  By  using 
the  Instrument's  peak  pick  routine  to  determine  absorption  values,  shifts  in 
the  absorption  baseline  caused  by  temperature  variations  in  the  room  were  not 
accounted  for.  When  the  absorbance  values  obtained  from  the  peak  pick 


routine  were  compared  with  those  obtained  from  manually  plotting  the  spectra, 
drawing  a  consistent  baseline  and  measuring  the  absorption,  it  was  clear  that 
a  different  approach  would  be  necessary  to  correctly  acquire  data  under 
instrument  control. 

To  properly  account  for  the  shift  in  baseline,  the  following  approach 
was  used.  For  each  molecular  species,  a  point  was  chosen  on  either  side  of 
the  analytical  wavelength,  essentially  the  same  points  one  would  choose  to 
draw  a  baseline  for  a  manual  determination  of  the  absorbance  from  a  plotted 
spectrum.  The  absorbance  values  at  these  two  points  were  determined  from  the 
instrument's  peak  pick  routine.  The  absorbance  value  at  the  analytical 
wavelength  was  also  determined,  as  before,  using  the  peak  pick  routine. 
However,  now  it  was  scaled  by  a  linear  interpolation  between  the  two  points 
on  either  side  (Reference  18).  This  procedure  was  implemented  by  using  a 
short  FORTRAN  program,  which  was  called  from  the  instrument's  MACRO 
software.  Tire  programs  are  listed  in  Appendix  B,  The  procedure  was  checked 
against  absorbance  values  obtained  by  manually  drawing  baselines  and 

measuring  absorbances.  The  results  were  essentially  identical.  Once  this 

latter  program  was  finished,  some  earlier  spectral  data  were  reprocessed  to 

correct  for  baseline  drift. 

This  procedure,  while  providing  the  correct  absorbance  values,  was 
operator-intensive.  Because  of  the  display- intensive  data  manipulation 
required  in  this  MACRO  program,  the  instrument  could  not  keep  up  with  the 
software  unless  distinct  software  pauses  were  inserted.  Each  pause  required 
the  operator  to  physically  press  the  RETURN  key  on  the  spectrometer  keyboard. 

To  alleviate  this  situation,  a  software  development  effort  was 

undertaken  with  the  goal  of  accomplishing  the  analysis  by  using  a  FORTRAN 
program.  The  basis  for  this  program  was  the  availability  of  an  existing 
FORTRAN  subroutine  called  FGETOF  written  by  Hicolet  which  allowed  FT-IR 
scratch  file  data  to  be  read  by  a  FORTRAN  program.  The  result  wae  a  program 
called  GET2  which  would  take  a  predetermined  wavelength  interval  on  one  aide 
of  an  analytical  peak  and  average  the  absorbance  values  over  that  interval. 
It  would  do  the  same  for  an  equal  interval  on  the  other  side  of  the 
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analytical  peak.  These  two  average  points  would  then  be  used  to 
baseline-correct  the  absorbance  at  the  analytical  peak. 

This  later  program  was  run  from  another  MACRO  (see  Appendix  C).  The 
program  GET2  was  used  to  correct  the  data  f4les  which  remained  uncorrected 
from  earlier  runs  (i.e.,  those  which  h*d  not  been  done  manually,  using  the 
earlier  MACRO  called  MPK) . 

3.  Oxygen  Concentration  Measurements. 

One  series  of  experiments  was  conducted  to  measure  the  permeation  of 
oxygen  into  the  chamber.  Oxygen  concentration  was  measured  by  an  indirect 
method,  using  a  dissolved  oxygen  (DO)  meter  (Yellow  Springs  Instruments, 
Model  58).  A  400-mL  beaker  containing  about  250  mL  of  water  was  placed  on  a 
magnetic  stirrer  in  the  chamber.  The  DO  meter  was  also  placed  in  the  chamber 
and  its  probe  (Yellow  Springs  Instruments,  Model  5720A)  was  supported  in  the 
water.  The  water  was  vigorously  stirred.  The  resulting  dissolved  oxygen 
reading  was  assumed  to  have  a  linear  proportionality  relationship  with  the 
concentration  of  oxygen  vapor  in  the  chamber.  The  DO  meter  readings  were 
recorded  continuously  on  a  strip-chart  recorder  (Esterline-Angus,  Model 
T171B).  The  chamber  mixing  fan  was  operated  continuously  during  these  oxygen 
concentration  measurements. 
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SECTION  III 


MODEL  DEVELOPMENT 


As  already  mentioned,  the  overall  reaction  of  gaseous  hydrazine  in  air  at 
room  temperature  up  to  160  °C  has  been  proposed  to  be: 


»2H4  +  °2  - ►  H2  +  2H2°  (1> 

In  addition  to  the  fact  that  chamber  surface-to-volume  ratios  affect  the  loss 
processes,  many  unexplained  experimental  results  in  hydrazine  studies  have 
also  been  obtained.  Most  notably,  some  results  show  hydrazine  approximately 
obeying  pseudo-first-order  kinetics*  others  clearly  do  not  (Reference  13). 


Studies  presented  in  this  report  have  been  geared  to  investigate  whether 
Equation  (1>  is  the  dominant  process,  and  if  not,  to  determine  what  kinetic 
processes,  chemical  or  physical,  are  occurring  and  to  successfully  model 
these  processes,  This  chapter  will  show  that  Equation  (1)  is  not  the  most 
important  gaseous  process,  and  will  develop  a  model  which  takas  into  account 
many  possible  processes;  including  gaseous  reactions,  surface-catalyzed 
reactions,  leakages,  gaseous  aerosol  formation,  surface  adsorption,  and 
permeation  through  chamber  walls, 


To  test  the  possibility  that  Equation  (1)  might  not  be  significant,  an 

experiment  was  conducted  in  which  the  Teflon  chamber  was  initially  filled 

with  dry  nitrogen  (fiL).  Hydrazine  was  then  vaporized  sad  flushed  into  the 

£  _ 

system  and  its  absorbance  was  monitored  with  time.  The  results  of  this  and 
other  experiments  showed  essentially  no  difference  in  the  hydrazine  loss  rate 
for  runs  in  !?2  and  in  air. 


Because  the  box  surrounding  the  Teflon  chamber  was  not  completely  oxygen 
(0,)-free,  oxygen  leakage  had  to  be  taken  into  account.  Before  the 

*4  .  ' 

hydrazine  experiments  (which  seem  to  show  no  sfta4or  oxygen  involvement)  could 
prove  to  be  conclusive,  a  model  describing  how  0^  slight  leak  into  the 
chamber  and  how  this  process  can  affect’  the  lass  rates  of  hydrazine  (if  it 
reacts  with  oxygen). was  developed. 


Two  possible  mechanisms  for  leakage  through  walls  are  effusion  and  perme¬ 
ation.  The  former  is  described  by  molecules  passing  through  holes  having 
cross-sectional  areas  that  are  small  when  compared  with  the  mean  free  path  of 
the  gas.  The  latter  is  described  by  molecules  jumping  from  one  active  site 
to  another  on  the  chamber  walls  or  by  simple  diffusion  caused  by  a  concentra¬ 
tion  gradient  (Reference  19).  Both  permeation  and  effusion  will  be 
incorporated  into  a  later  model.  Described  below  is  the  kinetic  model 
for  a  simple  effusion  process.  Although  results  shown  in  the  next  chapter 
indicate  permeation  is  probably  a  significant  process  in  Teflon  chambers,  the 
effusion  model  is  sufficient  to  account  for  gas  leakage  into  the  chamber. 

The  number  of  collisions,  v,  per  unit  area  per  unit  time  for  an  ideal  gas 

is: 

v  *  C\/RT/(2*M),  (2) 

in  which  C  is  the  molecular  concentration,  M  is  the  molecular  weight,  R  is 
the  gas  constant,  and  T  is  the  absolute  temperature.  The  number  of  molecules 
passing  through  an  aperture  of  area,  A,  in  time  At  is  CA\/Rt7'(2xH)! At .  If 
this  quantity  is  divided  by  NV,  in  which  N  is  Avogadro's  number  and  V  is  the 
chamber  volume,  then  the  resultant  quantity  would  be  the  change  in  moles  per 
unit  volume,  or  the  change  in  molar  concentration.  Tf  @  leaking  chamber  is 
initially  evacuated  of  a  gas  and  there  is  an  infinite  reservoir  of  that  gas 
outside  the  chamber  (Figure  6),  then  the  change  in  the  molar  concentration 
inside  the  chamber  is: 

ACm  =  (CffiQ  -  Cffl) ( A/V)\/RT/(2JtM) At ,  (3) 

in  which  C  is  C/tt,  the  chamber  molar  concentration,  and  <1  is  the 
m  mo 

reservoir  molar  concentration.  Taking  the  limit  as  At  approaches  aero  and 
integrating  yields: 

Cm  “  W1  ~  «sp[-(A/V)\/RT/(2aM)]t}.  <4> 

If  the  chamber  is  Initially  filled  with  the  gas  and  the  infinite  reservoir 
has  uone  of  the  gas,  C  can  be  shown  to  be: 
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Figure  6.  Diagram  Showing  Passage  of  Gas  Through  a  Hole  in  the 
Chamber  Wall. 
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Cm  =  Cmoexp [“<A/V) V  RT/<2«M) ] t  * 


(5) 


in  which  C^Q  is  now  the  initial  chamber  concentration.  Equations  (4)  and 

(5)  strictly  apply  if  the  mean  free  path  is  large  compared  to  the  cross 

section  of  the  aperture.  This  insures  that  all  molecules  passing  through  the 

holes  will  not  be  deflected  back  into  the  chamber  via  collisions.  At  25  °C 

—ft 

and  1  atm,  the  mean  free  path  of  0^  is  7  x  10”  meter,  which  is  much 
smaller  than  most  holes  which  might  exist  in  the  chamber  walls.  However, 
assuming  that  a  constant  fraction  of  molecules  passing  through  the  holes  are 
not  deflected  back  into  the  chamber,  the  concentration  change  can  still  be 
described  exponentially. 

The  differential  rate  equation  describing  the  process  in  Equation  (1), 
assuming  simple  second-order  kinetics,  is: 

d[N2H4]/dt  =  -k[N2H4] [02] ,  (6) 

in  which  k  is  the  second-order  rate  constant.  If  the  chamber  is  initially 
evacuated  of  02  and  hydrazine  is  promptly  injected  into  the  chamber,  then 
Equations  (4)  and  (6)  can  be  combined  to  yield  an  applicable  equation, 

d[N2H4J/dt  =  -k[N2H4][02]oCl  -  exp(-ket)],  (7) 

in  which  k  is  the  effusion  rate  constant  and  [0„]  is  the  reservoir 
e  2  0 

02  concentration.  Equation  (7)  assumes  the  initial  chamber  02 

concentration  is  zero  and  that  reaction  with  hydrazine  is  negligible  reletive 

to  02  effusion  into  the  chamber.  Since  [02]Q  >>  ^2H4^’  fckis 

latter  approximation  is  good  so  long  as  k  is  sufficiently  large. 

6 

Integrating  Equation  (7)  yields: 

[H2H4]  =  [N2H4]oexp(-kl02]o{t  -  [1  -  exp(-ket)]/ke>),  (8) 

in  which  [H2H43o  ia  the  initial  concentration  of  hydrazine.  Figure  7 

shows  plots  of  H  ]  vs  t  for  the  cases  k  =  k{0,l 

=  .05  hour”1,  and  k[0„]  =  .05  hour”1  and  k  =>  infinity;  and  a  plot 

Z  0  a  6 

of  [0J/[0J  vs  t  for  k  o  .05  hour”1.  If  Equation  (1)  describes 

2  2  0  £i 
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the  predominant  kinetic  process  for  hydrazine  in  the  Teflon  chambers,  then 
curve  1  of  Figure  7  should  depict  the  general  shape  if  the  chamber  is  ini¬ 
tially  filled  with  N.  and  k[0.]  is  comparable  to  or  larger  than  k  . 

Figure  8  shows  data  points  for  an  experiment  in  which  [02]  was 

monitored  indirectly  as  a  function  of  time  following  evacuation  of  02  from 

the  chamber.  The  0,  concentration  approximately  follows  curve  3  of  Figure 

z  -1 

7.  The  effusion  rate  constant  is  0.105  (.005)  hour  .  The  details  of  this 
experiment  are  given  in  the  Experimental  Section. 

Figure  9  shows  data  points  for  experiments  in  which  hydrazine  was 

monitored  in  the  chamber  filled  with  dry  air  and  in  the  chamber  initially 
filled  with  dry  N2>  As  indicated  previously,  there  is  little  difference  in 
the  plots.  These  two  experiments  were  conducted  in  the  same  S-mil  bag  just 
days  apart.  Hence,  physical  factors  which  may  have  a  gradual  affect  on 
changing  kinetic  behavior,  such  as  the  creation  of  more  holes,  were 
minimized.  In  addition,  the  fact  that  the  points  follow  a  general  decay 

pattern  and  not  a  sigmoidal  pattern  predicted  by  Equation  (8)  lends  further 

evidence  that  hydrazine  is  not  reacting  with  02«  The  decay  processes  must, 
therefore, be  one  or  more  of  the  following:  effusion,  permeation  through  the 
chamber  walls,  simple  physical  adsorption,  or  surface-catalyzed  reactions. 
Since  it  is  difficult  to  postulate  how  Teflon  surfaces  might  catalyze 
reactions,  one  or  more  of  the  other  processes  must  be  predominating.  Before 
the  idea  that  surface  interactions  are  important  can  be.  tested 
quantitatively,  a  comprehensive  model  which  fits  the  data  must  be  developed. 

The  following  kinetic  model  considers  all  likely  gaa  and  surface  inter¬ 
actions  which  might  occur  in  a  Teflon  chamber: 

ke 

6 - »•  effusion 

k 

G  +  8X  G-S! 

k-s 


2k 
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For  i  =  1  to  f; 


ki  . 

G  +  FjV  '  iG-Fj 
k_i 


For  i  =  1  to  p; 

kpi 

G  +  Pj  - reaction  products,  (9) 

in  which  G  is  the  gas  species  being  monitored,  is  the  first  polymeric 
film  pore  layer  (defined  below)  on  the  inside  surface  of  the  Teflon  chamber 
wall,  F^  is  the  Ith  interacting  surface,  other  than  Teflon,  is  the  Ith 
gaseous  species  which  reacts  irreversibly  in  a  simple  second-order  process 
with  G,  f  is  the  number  of  F-type  surfaces,  and  p  is  the  number  of  P-type  gas 
species.  The  corresponding  rate  constants  are  self-explanatory.  F^  can 
also  be  a  gas  species,  such  as  H20,  which  reversibly  complexes  with  G. 

Here  a  pore  layer  is  defined  as  the  distance,  perpendicular  to  the  plane  of 
the  wall,  that  a  molecule  must  travel  to  move  from  one  pore  site  to  an 
adjacent  site.  A  pore  layer  which  comprises  the  surface  of  the  film  is 
called  a  surface  pore  layer,  and  a  pore  layer  inside  the  film  is  called  a 
subsurface  pore  layer  (Figure  10).  There  are  two  surface  pore  layers,  one 
comprising  the  inside  surface  and  one  comprising  the  outside  surface  of  the 
film.  In  addition,  there  are  numerous  subsurface  pore  layers.  The  number  of 
pore  layers  is  that  number  comprising  the  thickness  of  the  walls. 
Delocalization  occurs  if  a  molecule  jumps  from  one  pore  site  to  another 
within  the  same  layer  (Figure  10). 

There  are  three  classifications  of  pores:  micropores  (less  than  15  X  in 
width),  mesopores  (15  to  1000  X  in  width),  and  macropores  (greater  than  1000 
X  in  width)  (Reference  20).  It  is  unlikely,  because  of  the  low 
concentrations  of  hydrazine  in  the  chamber,  that  more  than  one  molecule  will 
occupy  a  pore  site  at  any  given  time.  Hence,  the  kind  of  pore  which  will  occur  i 
Teflon  is  irrelevant  to  the  model.  Moreover,  whether  or  not  delocalization 
occurs  is  also  irrelevant  since  it  will  not  affect  the  model.  After  the  gas 
molecules  become  adsorbed  onto  the  surface  pore  layer,  permeation  and 
eventual  loss  can  occur  by  the  following  process: 
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Figure  10.  Diagram  Showing  How  Permeation  Can  Occur.  Molecules  Can 
Migrate  From  One  Layer  to  Another  (Arrow  1),  or  Migrate 
Within  the  Same  Layer  (Arrow  2) . 
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For  i  =  2  to  N; 


kD  t 

G“S^_2  +  Si  ^  ■  *  G— Si  +  Si_i 

kD 

k_g 

G-Sjj  — - — ►  loss  to  the  outside,  (10) 


in  which  S^  is  the  Ith  pore  layer,  k^  is  the  diffusion  rate  constant,  and 
N  is  the  number  of  pore  layers.  The  differential  rate  equations  describing 
the  processes  shown  in  Equations  (9)  and  (10)  are: 


f  P 

d[G]/dt  =  -(ke  +  ks[S]  +  I>itfi]  +  £kpi[Pi])[G]  +  k.gtG-Si]  + 

i=l  i=l 


f 

ZMte-Fii 

i=i 


d[G-Si]/dt  =  ks[S] [G]  -  k.alG-Sil  +  kD[S]([G-S2]  -  [G-Si] ) 

For  i  =  2  to  R-lj 

d[G-SiJ/dt  =  kuCSHtG-Si.x]  +  [G-Si+1]  -  2[G-Si]), 
d(G-SN]/dt  =  kDlSKtG-SH.!]  -  [G-SH])  -  k_s[G-SH] 

For  1  a  l  to  f; 

d[G-Fil/dt  =  kilFiHG]  -  k_i[G-FiJ,  (11) 

in  which  the  brackets  refer  to  the  concentrations,  and  [S]  is  the  concentra¬ 
tion  for  each  pore  layer.  If  [G]  becomes  comparable  to  [F^],  then  the 
pseudo-first-order  term,  kj(Fj][G],  must  be  replaced  by  the  second-order 
term,  ~  [G-Fj])[G],  in  which  (Fjlo  is  the  initial 

concentration  for  the  ltli  kind  of  surface.  All  surface  concentrations,  [S] , 
[Fj],  [G-S^J,  and  [G-F^],  are  in  the  same  units  as  the  gas  concentra¬ 
tions,  [Gj  and  [P^J,  and,  for  convenience,  are  defined  as  absorbance  units, 
which  is  proportional  to  the  gas  phase  concentration. 


29 


In  general,  the  series  of  differential  equations  above  cannot  be  analytically 
integrated.  However,  they  can  easily  be  numerically  solved  in  the  following 
manner.  If  the  variables,  x  and  t,  in  the  general  expression,  ax  =  g(x,t)At, 
cannot  be  separated,  then  each  infinitesimal  increment  in  x  can  be  calculated 
by  the  expression: 

x'  =  x  +  g(x,t)At.  (12) 

All  one  needs  to  know  from  the  physical  problem  are  the  initial  values  for  x 
and  t.  Each  increment  in  t  is  given  by: 

t’  =  t  +  At.  (13) 

A  computer  routine  can  be  written  to  calculate  each  incremental  value  for  x 
until  the  desired  value  of  t,  the  independent  variable,  is  reached.  The 
routine  requires  initial  conditions,  determined  by  the  physical  problem,  and 
a  value  for  At.  No  certain  value  for  At  is  required,  provided  it  is  small 
enough  to  insure  the  numeric  integration  converges.  As  it  becomes  smaller, 
more  computer  time  is  required  to  complete  the  integration. 

Although  this  technique  has  been  illustrated  for  just  one  differential 
equation,  it  can  in  theory  be  applied  to  an  infinite  number  of  dependent 
differential  equations.  Hence,  the  model  which  has  been  presented  can  be 
modified  and  expanded  to  any  degree  and  still  be  solved  numerically,  provided 
the  computer  can  do  the  integration  in  reasonable  time. 

This  model  can  be  fitted  to  actual  data  by  the  same  methods  used  to  fit 
simple  expressions  to  data.  The  only  modification  is  that  the  estimated 
values  for  the  dependent  variable  are  generated  from  a  numeric  Integration 
rather  than  a  simple  function.  The  parameters,  in  this  case  the  rate 
constants,  are  adjusted  to  give  the  best  estimates  in  the  same  manner.  The 
Harquardt  nonlinear  least  squares  approach  (run  on  a  Tektronix  computer, 
Model  4054)  was  used  to  fit  the  models  to  the  data.  The  program  used  was 
Statistics  Volume  4,  Tape  1.  This  program  allows  the  user  to  define  the 
function  and  set  initial  values  for  the  parameters.  Although  the  program  is 


in 


not  designed  for  the  user-defined  function  to  be  a  series  of  differential 
equations  requiring  numeric  integration,  the  program  is  able  to  fit  the 
differential  equations  to  the  data.  In  place  of  a  single  equation,  a 
user-written  subroutine  which  iteratively  performs  the  numeric  integration  is 
required.  The  subroutine  may  vary,  depending  on  how  the  model  is  used. 
Appendix  IV  lists  the  source  code  for  the  permeation  version  of  the  model  and 
explains  the  code. 

The  experimental  data  presented  in  this  report  include  monitoring  of 
temperature,  but  not  temperature  control.  At  the  beginning  of  this  study,  it 
was  thcught  that  temperature  fluctuations  affected  the  decay  of  hydrazine, 
hence  a  method  for  incorporating  temperature  fluctuations  into  the  model  was 
devised.  After  making  corrections  due  to  spectral  baseline  shifts  caused  by 
temperature  fluctuations,  it  was  found  that  temperature  fluctuations  had  no 
effect.  Nevertheless,  the  method  for  accounting  for  temperature  fluctuations 
is  outlined  below,  as  it  may  be  Important  for  some  processes  yet  to  be 
studied.  Because  the  equations  are  numerically  integrated,  it  is  an  easy 
task  to  include  temperature  variations  in  the  model. 

For  moderate  temperature  changes,  temperature  effects  can  be  incorporated 
into  the  rate  constants  by  the  Arrhenius  expression; 

k  .  A  exp(-Ea/(ET)J,  (14) 

in  which  A  is  the  pre-exponential  factor,  is  the  activation  energy,  and 
T  is  the  temperature.  In  order  for  the  differential  equations  to  be 
integrated  with  time,  the  temperature  must  be  expressed  as  a  function  of 
time.  The  general  expression  used  for  this  is: 

T  =  Tjil  ♦  f(t)J,  (15) 

in  which  Tj  is  the  temperature  where  f(t)  *  0  and  f(t)  is  some  function  of 
time,  t.  The  use  of  the  subscript  "l"  in  will  become  apparent  below. 
Since  it  Is  true  in  this  report  that  1  >>  f(t),  Equation  14  can  be  rearranged 
using  Equation  (15)  to  yield 
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(16) 


k  -  koexp[Eafvt)/(RT1)] , 


in  which  k  is  kQ  where  f(t)  =  0.  The  functional  form  of  f(t)  is  empiri¬ 
cal.  Generally,  though  not  always,  T  varies  sinusoidally,  being  a  maximum 
around  4  p.m.  and  a  minimum  around  7  a.m.  in  the  Environics  Laboratory. 

The  function,  f(t)  is  generated  one  of  two  ways.  One  way  is  to  fit  T  to 
an  expression  such  as  a  polynomial.  However,  T  can  sometimes  vary  in  some 
parts  of  the  kinetic  run,  but  remain  constant  over  several  hours  in  other 
parts.  In  these  cases  a  simple  empirical  equation  will  not  fit.  What  can  be 
used  instead  will  hereafter  be  called  "finite  linear  functions."  In  this 
approach,  T  is  assumed  to  vary  linearly  in  t  between  any  two  consecutive  data 
points.  This  approach  will  be  accurate  so  long  as  drastic  temperature 
fluctuations  do  not  occur  between  points.  The  procedure  for  applying  this 
approach  to  the  model  is  a3  follows.  Let  T^  and  t^  be  the  temperature 
and  time  for  the  Ith  data  point,  respectively.  Now  let  be  given  by  the 
following  general  linear  function: 

Ti  =  Tl^  +  mi t  i  +  (17). 

in  which  and  are  constants  associated  with  the  Ith  data  point. 
Since  temperature  varies  linearly  between  t^  and  ti+l'  m^,  b^ ,  and  T 
can  therefore  be  put  in  terms  of  T^,  t^,  tj+^,  ^i+i’  *  as 

follows: 


bi  -  Vh  - 1  -  -  V1 


”1  -  <T1+1  -  V'  'Vh.i  -  V1 


"  -  +  <h+i  -  V(t '  -  ti)- 


Equation  (18)  applies  in  the  range,  t^  i  t  f  t^  .  Therefore,  f(t),  in 


the  range  tj  £  t  £  ci+1,  is: 


Temperature  fluctuations  have  been  as  high  as  6.5  °C  during  the  course 
of  an  experiment.  Still  no  significant  effect  due  to  temperature  fluctua¬ 
tions  was  seen.  Hence,  no  processes  having  even  moderate  activation  energies 
are  occurring. 

Several  kinetic  reaction  models  which  do  lend  themselves  to  analytical 
integrations  have  been  derived  in  this  study.  Some  of  these  models  can  be 
found  in  many  advanced  books  on  chemical  kinetics.  Appendix  V  shows  the 
derivations  of  five  kinetic  models.  Some  of  them  have  been  used  to  fit  data 
from  this  study.  However,  they  did  not  fit  the  data  well. 

The  oxygen  permeation  data  was  modeled  using  an  adaptation  of  the 
permeation  version  of  the  model,  which  is  given  in  Appendix  VI. 

Some  experiments  were  done  in  which  Teflon  sheets  were  placed  inside  the 
chamber.  Appendix  VII  shows  an  adaptation  of  the  permeation  model  version 
which  accounts  for  permeable  sites  inside  the  chamber.  This  adaptation  has 
been  used  to  fit  data  from  the  runs  with  added  Teflon  sheets. 
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SECTION  IV 


RESULTS  AND  CONCLUSIONS 

A.  SPECIES  CONCENTRATION  MEASUREMENTS 

All  species  used  in  this  study  obeyed  Beer's  law  in  the  concentration 
ranges  used.  Kinetic  information,  therefore,  was  processed  by  using 
absorbance  versus  time  measurements. 

B.  THE  STABILITY  OF  INTERNAL  STANDARDS  IN  THE  TEFLON  CHAMBER 

Both  methane  (CH^)  and  sulfur  hexafluoride  (SF^)  were  used  to  test 
the  stability  of  inert  gases  in  the  Teflon  chamber.  Numerous  experimental 
runs  were  conducted  where  these  gases  were  introduced  into  the  chamber  under 
a  variety  of  experimental  conditions.  Methane  was  used  as  the  internal 
standard  for  most  of  the  hydrazine  runs  so  that  any  major  changes  in  chamber 
leak  rate  could  be  detected  and  repaired. 

i.  sr6 

A  total  of  7  different  runs  were  made  with  SF,  between  August  1985 

U 

and  September  1986.  These  Included  six  runs  in  dry  air  and  one  run  in  humid 
(50  percent  relative  humidity)  air.  The  concentration  of  SFfi  was 
approximately  equal  to  1.5  parts-per-million  by  volume  (ppm).  The  SF^  was 
injected  at  atmospheric  pressure  in  a  1.0-oL  gas  syringe  as  a  0.5  mL  sample. 
All  of  the  runs  in  both  dry  and  humid  air  were  made  in  combination  with 
approximately  20  ppm  CH^. 

The  decay  of  SF,  was  erratic,  A  typical  decay  curve  (a  plot  of 
o 

SF,  absorbance  versus  time)  is  shown  in  Figure  11.  Further  investigation 
o 

of  this  behavior  revealed  that  the  majority  of  the  peaks  and  valleys  in  the 
decay  curve  decreased  considerably  when  the  integrated  band  intensity  rather 
than  the  absorbance  was  plotted  versus  time  (see  Figure  11).  The  revised 
data  processing  MACRO  program  is  listed  as  a  part  of  Appendix  B.  The  erratic 
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Figure  11.  Comparison  of  SP  Plots  by  Integration  and  Absorbance. 


behavior  in  the  plots  is  understandable  in  view  of  the  fact  that  the 
absorption 

hydrazine.  Plotting  SF, 


absorption  coefficient  for  SF,  is  about  100  times  greater  than  that  of 

o 


absorbance  versus  time  has  the  effect  of 
amplifying  small  variations  by  a  factor  of  100  when  compared  with  a  similar 
hydrazine  decay  curve. 


The  SF,  runs  were  analyzed  by  plotting  the  integrated  absorbance  of  the 
6  -1 

SFg  band  between  900  and  1000  cm  versus  time  (normalized  to  a  starting 
value  of  one  for  comparison).  They  are  shown  in  Figure  12.  The  decay 
curves  fall  roughly  into  two  groups.  The  upper  four  curves  have  an  average 
half-life  of  about  60  hours.  The  lower  three  curves  have  an  average 
half-life  of  about  40  hours. 


An  examination  of  the  plot  legend  shows  that  the  upper  four  curves  were 
obtained  in  experiments  from  October  1985  through  February  1986.  The  lower 
three  curves  were  from  experiments  from  July  1986  through  September  1986. 
The  33  percent  decrease  in  half-life  is  probably  due  to  leaks  from  small 
holes  created  by  mechanical  stress  on  the  Teflon  bag  as  a  routine  part  of 
other  experiments  conducted  from  March  1986  to  June  1986.  Repeated 
experiments,  especially  those  where  the  bag  was  opened,  often  caused  small 
holes  along  seams  and  at  corners. 


2.  CH. 


Methane  was  used  as  an  internal  standard  or  an  inert  test  gas  in  22 
separate  experimental  runs  between  10  December  1985  and  19  September  1986. 
Five  of  these  runs  were  conducted  in  humid  air  and  the  remainder  in  dry  air. 
Six  of  the  runs  contained  SF,,  three  contained  ammonia  (NH«),  and  one 

6  J 

contained  additional  Teflon  sheets. 

In  the  runs  made  with  SF&,  when  normalized  to  a  beginning 
absorbance  of  one,  the  CH^  decay  followed  the  same  pattern  as  the  SF^. 
These  curves  are  plotted  in  Figure  13.  There  were  roughly  two  groups;  those 
conducted  early  in  the  year  had  a  half-life  of  roughly  60  hours  and  those 
conducted  from  July  1986  to  September  1986  had  a  half-life  of  35  to  45  hours. 
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Figure  i2.  Plots  of  Stv  Concentration  Versus  Time  in  the  Teflon  Chamber 
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Figure  13.  Plots  of  CH^  Concentration  Versus  Time  in  the  Teflon  Chamber 
(with  SF&  Prese:  t) . 
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The  single  run  conducted  under  humid  conditions  showed  a  somewhat  slower 
decay  rate  than  the  runs  conducted  under  dry  conditions.  This  pattern 

appears  throughout  the  CH^  runs  and  demonstrates  that  either  there  is  a 
site  competition  between  CH^  and  water  or  that  the  humid  air  is  wetting  the 
Teflon  film,  making  it  less  receptive  to  CH^  adsorption. 

In  experiments  conducted  with  NH^,  the  run  conducted  under  humid 
conditions  decayed  somewhat  more  slowly  than  the  two  conducted  under  dry 
conditions  (half-life  42  hours  versus  35  hours).  Even  though  these  runs  were 
all  conducted  early  in  the  year,  they  did  not  follow  the  trend  of  slower 
decay  rates  observed  in  the  SF,/CH.  system.  The  reason  for  this 

0  4 

difference  is  not  clear.  These  three  runs  are  plotted  in  Figure  14. 

In  the  experiments  where  hydrazine  was  run  in  combination  with  CH^, 

the  decay  rates  generally  followed  the  same  pattern  observed  in  the 

SF,/CH.  system.  Decay  rates  were  generally  slower  in  the  runs  conducted 
6  4 

at  the  first  part  of  the  year  and  somewhat  faster  in  runs  conducted  some 

months  later.  Again,  the  runs  conducted  under  humid  conditions  decayed 

somewhat  more  slowly  than  any  of  the  others.  These  runs  are  shown  in  Figure 

15.  CH,  runs  conducted  with  added  Teflon  sheets  and  added  mirror  surfaces 
4 

were  not  distinguishable  from  others. 

One  other  observation  should  be  noted.  In  runs  conducted  under  humid 
conditions  (with  the  exception  of  the  NH^/CH^  run),  the  CH^  concentra¬ 
tion  increased  for  the  first  4  to  6  hours  of  the  run.  It  would  then  show  a 
smooth  gradual  decrease.  This  behavior  is  shown  in  Figure  16. 

3.  Other  Test  Gases 

Other  runs  were  conducted  with  NH^,  methanol,  dichlorodifluoro- 
methane  (Freon-12),  atmospheric  0^,  and  water  vapor  to  test  ideas 
concerning  properties  .of  the  Teflon  film. 
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Figure  14.  Plots  of  CH4  Concentration  Versus  Time  in  the  Teflon  Chamber 
(with  NH,  Present) . 
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Figure  15.  Plots  of  CH^  Concentration  Versus  Time  in  the  Teflon  Chamber 
(with  Hydrazine  Present) . 
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Figure  16.  Plots  of  CH^  Concentration  Versus  Time  in  the  Teflon  Chamber 
(under  High  Humidity  Conditions) . 
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a.  NH3 

Since  hydrazine  may  be  considered  as  being  derived  from  NH^  by 
replacement  of  one  hydrogen  atom  with  an  -NH2  group  (Reference  21),  it  was 
felt  that  NH^  would  exhibit  behavior  similar  to  hydrazine  in  the  Teflon 
chamber.  Three  different  experiments  were  run  with  NH^,  two  in  dry  air  and 
one  in  humid  air.  In  each  case,  CH^  was  used  as  an  internal  standard. 

These  runs  were  normalized  to  an  initial  integrated  absorbance  of 
one  and  plotted  in  Figure  17.  The  NH^  decay  in  humid  air  was  a  little 
slower  than  in  dry  air.  The  estimated  half-life  in  humid  air  was  about  41 
hours,  while  the  average  half-life  in  dry  air  was  about  36  hours. 

In  general,  NH^  behaved  in  the  same  way  as  the  inert  tracer 
gases  had  behaved.  Despite  its  chemical  similarity  with  hydrazine,  its 
behavior  in  the  Teflon  chamber  failed  to  indicate  any  hydrazine-like  decay. 

b.  Methanol 

Since  NH3  behaved  essentially  the  same  as  CH^  or  SFg  in  the 
Teflon  chamber,  methanol  was  tried  in  one  run  to  see  if  it  would  behave  more 
like  hydrazine.  The  run  was  conducted  in  dry  air,  and  methanol,  too,  behaved 
like  the  inert  gases.  The  run  is  plotted  is  Figure  18. 

c.  Freon-12 

Technical  data  on  the  permeation  rates  of  different  gases  through 
FEP  Teflon  film  were  obtained  from  DuPont  (Reference  22).  The  permeation 
rate  of  Freon-12  was  from  30  to  100  times  greater  than  most  other  gases.  One 
experiment  was  conducted  with  Freon-12  in  dry  air  to  determine  its  behavior 
in  the  Teflon  chamber.  The  expectation  was  that  the  Freon-12  would  rapidly 
permeate  out  of  the  chamber  and  have  a  relatively  short  half-life.  However, 
it,  too,  behaved  much  the  same  as  Cli^  or  SF^  (see  Figure  19). 
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Figure  17.  Plots  of  NH^  Concentration  Versus  Tine  in  the  Teflon  Chamber. 


The  behavior  of  Freon-12  in  the  chamber  is  difficult  to  explain 
unless  one  invokes  the  idea  that  some  gases  in  the  chamber  escaped 
predominantly  from  undetected  macroscopic  holes  in  the  chamber  walls  or 
seams.  This  possibility  cannot  be  totally  ruled  out,  since  small  holes  are 
very  difficult  to  locate  in  the  transparent  Teflon  film  surface,  especially 
at  corners  and  seams, 

d.  Water  Vapor 

Water  permeation  studies  were  done  in  which  the  chamber  was 
filled  with  dry  air  and  the  increase  in  water  vapor  was  monitored  with  the 
FT-IR  spectrometer.  Figure  20  shows  a  first-order  fit  to  a  typical  data 
set.  The  half-lives  estimated  from  the  various  data  sets  varied  somewhat. 
The  data  in  Figure  20  showed  a  half-life  of  19  hours.  Most  half-lives  were 
smaller  than  this.  Moreover,  some  data  sets  showed  substantial  scatter  after 
15  to  20  hours.  This  may  have  been  due  either  to  temperature  fluctuations  or 
changes  in  the  absolute  humidity  in  the  laboratory. 

Water  vapor  has  a  somewhat  greater  ability  to  go  through  Teflon 
film  than  most  of  the  other  gases  studied.  That  water  may  have  a  substantial 
affinity  for  Teflon  film  is  an  important  consideration  in  light  of  the 
hydrazine  studies  presented  later  in  this  report. 

e.  Atmospheric  0^ 

Oxygen  permeation  studies  were  done  in  which  the  chamber  was 
filled  with  dry  N2  and  the  increase  in  02  was  monitored  with  the  FT-IR 
spectrometer,  as  described  earlier.  The  estimated  half-life  for  this  process 
was  7  hours.  This  process  was  modeled  using  the  simple  effusion  model  in 
Section  III,  Model  Development.  Figure  21  shows  a  fit  of  the  permeation 
version  adapted  tor  permeation  into  a  chamber  (see  Appendix  F).  The  fitted 
parameters  with  errors  are:  kg [ SJ  »  ,7(.2),  k^CSl  «  2('>),  and  = 
.8(.4)  (all  parameters  in  hour”1).  This  model  aptly  explains  the  slight 
sigmoidal  shape  of  the  data.  Oxygen  behaves  much  like  hydraaine  toward 
Teflon  film,  despite  the  obvious  structural  and  chemical  differences. 
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Figure  20.  First-Order  Fit  Co  Water  Permeation  Data. 
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f.  Discussion 


There  are  two  modes  by  which  a  gas  can  enter  or  leave  the 
chamber:  permeation  and  effusion.  Experiments  with  inert  gases  conducted 
over  a  period  of  months  seem  to  indicate  that  effusion  is  occurring  to  a 
larger  extent.  However,  if  effusion  was  the  only  mode  of  loss  from  the 
chamber,  there  would  be  a  correlation  between  the  square  root  of  the  mass  of 
a  chemical  species  and  its  loss  rate.  No  such  correlation  was  observed  in 
this  study.  Therefore,  permeation  must  also  be  an  important  loss  mode.  The 
degree  of  importance  of  permeation  may  be  somewhat  dependent  on  the  chemical 
nature  of  the  species.  However,  there  is  no  apparent  correlation  between  the 
structural  and  chemical  nature  of  a  species  and  its  loss  rate.  There  may  be 
some  less  obvious  structural  or  chemical  feature,  either  inherent  or  induced, 
common  to  hydrazine,  water,  and  oxygen,  which  accounts  for  the  enhanced 
permeation  rates  of  these  species. 

The  fact  that  the  0^  permeation  data  can  be  fitted  very  well  by 
the  permeation  version  of  the  model  may  seem  to  lend  support  to  the 
hypothesis  that  permeation  is  important  for  some  species  in  Teflon  film 
chambers.  However,  a  model  version  in  which  simple  effusion  is  coupled  with 
surface  adsorption  onto  nonpermeable  sites  inside  the  chamber  fits  the  0^ 
data  equally  well.  Experimental  results  from  all  inert  species  studied 
indicate  both  effusion  and  permeation  are  taking  place,  but  the  model  cannot 
be  used  to  determine  the  relative  importance. 

C.  HYDRAZINE  EXPERIMENTS 

Nearly  40  r  jparate  runs  were  conducted  with  hydrazine  in  the  5-mil  Teflon 
chamber  beginning  in  October  1985.  These  included  runs  in  dry  air,  in  humid 
air,  in  dry  N.,,  and  in  humid  Nj.  There  were  also  runs  made  with  added 
Teflon  sheets  and  with  added  mirror  surfaces  in  the  chamber. 

Various  abridged  versions  of  the  model  have  been  used  to  fit  single  runs 
and  sequential  runs.  Versions  used  for  single  runs  included  the  permeation 
version  and  the  coupled  effusion-adsorption  version.  These  simple  versions 
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fit  data  sets  much  better  than  the  first-order  equation,  but  do  not  account 
for  conditioning.  A  simple  modification  of  the  effusion-adsorption  version, 
however,  can  account  for  conditioning.  Changing  the  adsorption  process  from 
first-  to  second-order  is  sufficient  to  handle  the  conditioning  which  is  seen 
in  sequential  data  sets. 

1.  Hydrazine  in  Dry  N 

A  total  of  five  different  experiments  were  conducted.  The  runs  were 
normalized  to  a  starting  absorbance  of  one  for  the  hydrazine  analytical  band 
at  958  cm  1  and  plotted  on  one  graph  in  Figure  22.  Four  of  the  runs  fell 
on  nearly  the  same  decay  curve,  while  one  (10-12  Jan  86)  is 
considerably  slower.  A  simple  graphical  determination  of  the  half-lives  for 
these  runs  gave  the  results  shown  in  Table  1. 


TABLE  1.  HYDRAZINE  DECAY  IN  DRY  NITROGEN 


Date  of  Run 


Half-Life, (hours) 


6-8  Dec  85  8.0 
11-12  Dec  85  8.1 
10-12  Jan  86  16.8 
19-21  Feb  86  6.9 
8-9  Sep  86  9.5 


In  examining  this  data,  it  was  found  that  the  10-12  Jan  86  run  had 
only  one  observable  difference  from  the  other  four  runs.  It  was  the  last  in 
a  series  of  four  experimental  runs  involving  hydrazine  decay  which  were 
conducted  essentially  back-to-back;  that  is,  there  was  never  more  than  a 
weekend  between  them.  Even  though  each  run  was  preceded  by  a  thorough  purge, 
there  is  apparently  enough  surface  site  attraction  for  hydrazine  to  maintain 
surface  conditioning  for  at  least  several  days.  Thus,  the  chamber  was 
conditioned  by  the  earlier  runs  and,  as  a  result,  the  decay  rate  was  only 
about  half  that  of  the  other  runs.  In  fact,  model  fits  shown  later  indicate 
that  deconditioning  is  a  relatively  slow  process  in  dry  conditions. 

Three  of  the  other  runs  had  been  conducted  after  lengthy  periods  of 
no  chamber  activity:  5  weeks  for  the  6-8  Dec  85  run,  2  weeks  for  the  19-21 
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Feb  86  run,  and  5  weeks  for  the  8-9  Sep  86  run.  However,  the  11-12  Dec  85 
run  was  conducted  only  3  days  after  the  6-8  Dec  85  run.  There  was,  however, 
a  dry  SF6/CH4  run  conducted  between  6-8  Dec  85  and  11-23  Dec  85.  This 
intervening  run  apparently  affected  the  chamber  deconditioning  process,  but 
it  is  not  clear  by  what  mechanism. 

Figure  23  shows  the  data  for  the  8-9  Sep  86  run  with  model  fits.  The 
permeation  version  fits  the  data  much  better  than  a  simple  first-order  fit. 
The  fitted  parameters  are:  kg[S]  =  .17(.01),  k_g  =  .58(.07),  and  kptS] 
=  «9(.2)  (all  parameters  are  in  hour  ^).  The  number  of  pore  layers  was  set 
to  two  and  the  time  step  was  .01  hour.  This  was  generally  the  case  for  most 
of  the  permeation  version  fits. 

2.  Hydrazine  in  Dry  Air 

Beginning  in  October  1985  and  continuing  through  December  1986,  a 
total  of  18  different  experiments  of  hydrazine  in  dry  air  were  conducted. 
Fifteen  of  these  runs  were  normalized  to  a  starting  absorbance  of  one  for 
comparison  with  each  other  (the  remaining  three  were  not  successfully 
stored).  These  runs  are  shown  in  Figures  24,  25,  and  26.  The  half-life  for 
each  run  was  again  estimated  by  a  simple  graphical  procedure.  The  results 
are  tabulated  in  Table  2. 

The  first  group  of  data  is  a  summary  of  experimental  runs  conducted 

between  October  1985  and  February  1986.  These  runs  are  plotted  in  Figure 

24.  The  28-29  Oct  85  run  was  the  first  hydrazine  run  in  the  5-mil  Teflon 

bag.  However,  two  characterization  runs  had  been  conducted  earlier,  one  on 

16  Oct  and  the  other  on  21  Oct.  These  runs  were  made  With  CH.  and  SF,  to 

4  6 

evaluate  the  stability  of  these  inert  gases  in  the  chamber.  Apparently  these 
two  runs  affect  the  chamber  conditioning  (giving  a  half-life  about  30  percent 
greater  than  later  runs  in  a  well-flushed  chamber),  but  the  exact  mechanism 
is  not  yet  established.  The  next  run,  conducted  from  19-20  Dec  85  showed  a 
half-life  of  only  7.8  hours.  This  decrease  is  explained  by  the  fact  that  the 
Teflon  chamber  sat  for  7  days  before  the  run  was  started,  giving  the  chamber 
time  to  partially  decondition.  The  third  run  in  the  table,  23-24  Dec  85,  was 
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Figure  24.  Plots  of  Hydr2i.n0  Concentration  Versus  Tima  in  the  Teflon  Chamber 
(In  Dry  Air). 
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TABLE  2.  HYDRAZINE  DECAY  IN  DRY  AIR 


Date  of  Run 


Half-Life  (hours) 


-  Group  1 


28-29 

Oct 

85 

10.0 

19-20 

Dec 

85 

7.8 

23-24 

Dec 

85 

10.8 

3-5 

Feb 

86 

6.6 

21-22 

Feb 

86 

6.5 

Group  2 


14-15  Apr  86 

9.6 

15-16  Apr  86 

12.9 

16-17  Apr  86 

14.0 

17-18  Apr  86 

14,6 

- Group  3  — - 


3-5  Nov  86 

6.7 

12-14  Nov  86 

8.6 

24-26  Nov  86 

10.4 

1-3  Dec  86 

10.3 

8-9  Dec  86 

10.3 

9-11  Dec  86 

8.7 

conducted  2  days  after  the  second  run.  This  did  not  allow  the  chamber  walls 
to  desorb  the  hydrazine  completely  and  the  half-life  increased  to  10.8  hours. 
Prior  to  the  fourth  run,  8  sheets  of  2-ail  Teflon  were  hung  from  the  interior 
celling  of  the  chamber,  effectively  doubling  the  interior  surface  area. 
However,  the  chamber  had  also  sat  for  2  weeks  before  the  run  was  made.  The 
combination  of  added  Teflon  sheets  and  a  2-week  rest  period  combined  to  give 
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a  very  short  6.6  hour  half-life.  The  final  run  in  the  first  group  was  made 
after  eight  Pyrex  mirrors  of  various  sizes  were  placed  inside  the  chamber. 
The  mirrors  were  of  the  same  basic  construction  as  those  used  in  the  chamber 
optical  system.  They  were  of  Pyrex  glass  with  a  silicon  dioxide  overcoated 
aluminum  surface.  The  run  was  made  to  test  the  effect  of  added  mirrors  on 
the  decay  of  hydrazine.  This  run  was  preceded  by  a  run  with  dry  nitrogen; 
therefore,  there  had  not  been  enough  time  for  the  chamber  surfaces  to  desorb 
hydrazine  completely.  The  added  mirrors  appear  to  have  accelerated  the  decay 
process,  since  the  half-life  was  the  lowest  of  any  run,  6.5  hours. 

The  third  run  in  Group  1  was  modeled  using  an  adaptation  of  the 
permeation  version  of  the  model.  This  adaptation,  derived  in  Appendix  G, 
includes  equations  for  a  permeable  surface  within  the  chamber.  Since  2-mil 
sheets  were  used  in  the  chamber  made  of  5-mil  sheets,  the  number  of  inside 
surface  pore  layers  was  set  to  two  and  the  number  of  chamber  wall  pore  layers 
was  set  to  five.  Table  3  shows  results  of  the  model  fit  to  this  data  set, 
along  with  three  more  fits  of  data  sets  in  dry  conditions.  Since  long-term 
conditioning  is  not  taken  into  account  in  the  permeation  version  of  the 

model,  the  fitted  parameters  shown  in  Table  3  are  for  rough  comparisons 
only.  The  value  of  kg(S]  for  the  run  with  the  added  Bheets  is  slightly 

higher  than  for  the  runs  with  nothing  added,  reflecting  the  idea  that  the 

2-mil  sheets  were  unconditioned.  The  parametric  values  for  the  fit  of  the 
28-29  Oct  85  run  are  somewhat  different  than  for  the  other  two  fits.  Since 
the  run  was  the  first  run  in  the  5-mil  bag,  this  discrepancy  might  be 
expected.  However,  because  the  permeation  version  is  not  an  accurate  model 
(it  does  not  take  into  account  conditioning),  differences  in  the  values  for 
the  fitted  parameters  cannot  be  confidently  interpreted. 

The  value  of  k  [S]  for  the  fit  of  the  run  with  the  added  mirrors  is 

s 

larger  than  for  the  other  fits.  This  lends  support  to  the  idea  that  the 

hydrazine  is  interacting  with  the  mirrors,  either  by  surface  adsorption  or 
surface-catalyzed  oxidation.  A  more  accurate  version  of  the  model  would  take 
these  interactions  into  account  by  introducing  an  additional  adsorption 
process,  but  would  also  require  additional  parameters.  As  it  is,  the  simple 
permeation  version  fits  the  data  quite  well.  Using  more  sophisticated 


TABLE  3.  RESULTS  OF  THE  PERMEATION  VERSION  MODEL  FITS  FOR  REPRESENTATIVE 
DATA  SETS.* 


Date  of  Run 

Conditions 

ka[S] 

k— 3 

kD[S] 

SOS 

28-29  Oct  85 

dry  air 

.213( .007) 

1.96( .09) 

5 . 5  ( .  2  ) 

3 

x  10“3 

3-5  Feb  86 

dry  air; 
2-mil  sheets 

.24( .03) 

3.2( .4) 

28(13) 

1 

I 

o 

H 

X 

21-22  Feb  86 

dry  N25 

8  mirrors 

•31( .02) 

2.1(.2) 

5.7( .4) 

9 

x  10-*4 

14-15  Apr  86 

dry  air 

.18( .03) 

3.3C.8) 

32(10) 

5 

x  10"4 

*The  model  was  adapted  as  described  in  the  text  for  the  data  set  using 
the  2-mil  Teflon  sheets.  The  number  of  pore  -ayera  for  the  chamber  film 
for  all  runs  was  set  to  5.  All  parameters  are  in  hour*1.  SOS  is  the 
residual  sum  of  acuares. 

versions  of  the  model  requires  fitting  multiple  data  sets  simultaneously 
rather  than  just  single  data  sets.  Fits  of  multiple-data  sets  have  been  done 
for  sequential  runs  and  the  results  are  given  later  in  this  section. 

The  second  group  of  data  in  Table  2  summarizes  four  sequential  runs  made 
in  April  1986.  These  runs  are  shown  in  Figure  25*  These  runs  were  made  to 
teat  the  hypothesis  that  there  was  a  conditioning  effect  when  insufficient 
time  was  al lowed  between  hydrazine  decay  runs.  The  first  run  in  this  series 
was  conducted  after  the  chamber  had  not  been  in  use  for  almost  5  weeks. 
According  to  previous  experience,  the  half-life  should  have  been  around  7 
hours.  However,  it  was  measured  at  9.6  hours.  The  reason  for  this  behavior 
is  not  known.  Each  succeeding  run  in  this  series  was  conducted  by  adding 
hydrazine  to  the  previous  run  with  no  additional  flushing  of  the  chamber.  As 
expected,  the  half-life  of  runs  1,  2,  3,  and  A  increased  from  9.6  to  14.6 
hours.  This  series  of  runs  demonstrated  that  there  is  a  definite 
conditioning  process  which  occurs  in  this  chamber. 

The  second  group  in  Table  2  was  modeled  compositely,  using  the  effusion- 
adsorption  version  modified  by  allowing  adsorption  to  he  second-order.  The 
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parameters  for  this  version  are:  k^  (adsorption  constant),  k  ^ 
(desorption  constant),  (initial  concentration  of  nonpermeable 
surface  sites),  and  Icg  (first-order  effusion  constant).  This  version  has 
been  used  for  composite  fits  of  sequential  runs  because  no  reasonable 
adaptation  of  the  permeation  version  could  fit  a  set  of  sequential  runs  very 
well.  This  is  not  to  say  that  permeation  does  not  occur,  at  least  to  some 
extent.  In  fact,  the  effusion  constant  serves  as  a  catch-all  term  for  any 
first-order  process,  which  includes  steady-state  permeaticn.  Figure  27  shows 
a  normalized  composite  fit  to  the  first  three  runs  of  the  sequential  data 
set.  Initial  absorbance  readings  for  the  second  and  third  runs  were  adjusted 
to  minimize  the  residual  sum  of  squares.  This  was  done  because  of  the  large 
systematic  errors  in  the  first  data  points  of  the  second  and  third  rums. 
Table  4  shows  results  of  the  composite  fits  for  the  second  and  third  groups 
of  hydrazine  decay  in  dry  air  and  two  composite  fits  for  runs  in  humid  air. 
Since  the  effusion-adsorption  version  adapted  for  describing  the  conditioning 
process  is  likely  a  more  accurate  description  of  the  truth  than  the  simple 
permeation  version,  a  meaningful  interpretation  of  the  fitted  parameters  is 
possible.  Interpretation  of  these  parameters  is  given  later  in  the  General 
Discussion  Section. 

The  third  group  of  data  in  Table  2  summarizes  six  sequential  runs  made 
between  3  Nov  and  11  Dec  86  (Figure  25).  These  runs  were  designed  to  further 
test  the  model.  The  results  of  the  composite  fit  are  in  Table  4.  The  first 
run  was  begun  after  a  2-week  purge  of  the  chamber  with  pure  air.  Every  30 
minutes,  the  mixing  fan  was  cycled  on  for  2  minutes  and  the  annular  space  in 
the  Plexiglas  box  was  continuously  purged  with  dry  air.  The  expectation  was 
that  the  initial  run  would  be  essentially  the  same  as  a  run  in  *>  "new"  Teflon 
chamber. 

This  turned  out  to  be  the  case;  the  half-life  was  6.7  hours.  Runs  2-5 
were  conducted  by  adding  additional  hydrazine  to  the  chamber.  There  was  no 
purging  between  runs,  but  the  chamber  mixing  fan  was  left  on  its  automatic 
cycle.  The  second  run  saw  an  increase  in  half-life  to  8.6  hours  and  the 
third  run  to  10.4  hours.  The  adsorption  sites  had  become  nearly  saturated 
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Figure  25.  Plots  of  Hydrazine  Concentration  Versus  Time  in  the  Teflon  Chamber 
(Sequential  Runs  in  Dry  Air) . 
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Figure  27.  Effusion-Second-Order  Adsorption  Model  Fit  to  a  Sequential 
Set  of  Hydrazine  Decay  Runs  in  Dry  Air. 
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after  the  third  run,  causing  the  half-lives  for  Runs  4  and  5  to  be  nearly  the 
same. 

Run  6  was  conducted  by  adding  hydrazine  to  Run  5  at  the  24-hour  point. 
The  resulting  decrease  in  half-life  is  somewhat  misleading.  A  better 
understanding  of  the  situation  in  this  run  is  obtained  by  looking  at  the 
shape  of  the  decay  curve  in  Figure  26.  The  model  predicted  that  the  decay 
curve  for  this  final  run  would  lie  just  a  little  above  that  of  Run  5. 
Instead,  Run  6  began  to  decay  more  rapidly  than  any  of  the  previous  runs  but 
then  ended  up  with  slower  decay  beyond  about  the  30-hour  point.  The 
explanation  for  this  decay  pattern  is  not  known. 

There  was  no  readily  observable  difference  between  the  runs  conducted  in 
dry  nitrogen  and  the  first  group  of  those  conducted  in  dry  air.  Succeeding 
groups  in  dry  air  showed  decay  half-lives  which  increased  to  an  essentially 
constant  level.  No  such  experiments  were  conducted  with  dry  nitrogen  because 
of  the  large  quantities  required  for  continuous  purging.  If  a  difference 
between  oxygen  and  nitrogen  decay  rate  is  to  be  observed  in  this  particular 
chamber,  it  would  probably  require  that  long-term  successive-addition 
hydrazine  decay  runs  be  made  in  a  dry  nitrogen  atmosphere. 

3.  Hydrazine  in  Humid  Nitrogen 

Only  a  single  experimental  run  was  conducted  under  these 
conditions.  The  decay  curve  is  shown  in  Figure  28.  Several  observations 
were  noted  in  this  system.  First,  the  decay  curve  is  not  pseudo-first- 
order.  (Figure  29  shov-r  results  of  a  first-order  fit  and  the  permeation 
version  fit).  Second,  the  decay  rate  is  much  more  rapid,  giving  an 
approximate  half-1' C-  of  about  3.3  hours.  Third,  there  is  a  fairly 
significant  amount  of  NH^  produced,  about  4  ppm. 

The  explanation  for  the  dramatic  differences  in  the  decay  curves  for 
this  system  versus  '•.lie  hydrazine-dry  nitrogen  system  appears  to  lie  in  the 
ability  of  hydrazine  to  readily  adsorb  onto  wetted  Teflon  film,  just  as  in 
the  case  of  humid  air.  In  all  likelihood,  all  of  the  "active  sites"  (both 
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Figure  28.  Plot  of  Hydrazine  Concentration  Versus  Tine  in  the  Teflon  Chamber 
(In  Humid  Nitrogen) . 
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surface  sites  and  microvoids  in  the  film  itself)  on  the  Teflon  chamber 
surface  are  occupied  by  water  when  the  hydrazine  is  added.  As  the  hydrazine 
molecules  strike  the  surface  of  the  chamber,  they  encounter  these  adsorbed 
water  molecules  or  microdroplets  and  stick.  The  presence  of  the  adsorbed 
water  apparently  exerts  a  catalytic  effect  on  the  trapped  hydrazine,  causing 
the  formation  of  small  quantities  of  NH^. 

4.  Hydrazine  in  Humid  Air 

A  total  of  12  runs  were  completed  in  humid  air.  The  first  two  were 
conducted  in  sequences  of  runs  involving  several  experimental  conditions. 
The  next  four  were  run  as  a  series  with  no  purging  of  the  chamber  between 
them.  The  last  six  were  conducted  after  6  weeks  of  purging  the  chamber  with 
humid  air  (with  the  mixing  fan  on  a  30-minutes-off ,  2-oinutes-on  cycle)  and 
the  Plexiglas  box  with  humid  air.  This  was  designed  to  give  conditions 
approximating  those  which  would  be  found  in  a  new  Teflon  bag  which  had  been 
conditioned  with  water  vapor.  The  first  two  runs  in  the  sequential  data  sets 
were  fitted  using  the  effusion-adsorption  version  adapted  for  conditioning 
and  results  are  given  in  Table  4.  The  dates  and  half-lives  of  these  runs  are 
shown  in  Table  S.  Plots  of  typical  data  sets  are  shown  in  Figures  30  and  31 
(decay  curves  are  normalized  to  unit  absorbance). 

The  30  Oct  85  run  was  the  first  one  in  the  5-mil  Teflon  chamber  involving 
high  humidity.  In  this  case,  room  air,  measured  at  74  percent  relative 
humidity,  was  used.  However,  this  run  was  preceded  by  a  hydrazine  run  in  dry 
nitrogen  on  28  Oct  and  a  hydrazine-conditioning  run  (i«e.,  an  injection  of 
four  times  the  regular  amount  of  hydrazine)  on  29  Oct.  Apparently  the  high- 
humidity  room  air  negated  the  expected  conditioning  effects  of  the  two 
preceding  runs,  resulting  in  the  very  short  half-life  observed  in  the  30  Oct 
experiment.  Model  fits  on  humid  runs  show  that  desorption  is  greatly 
enhanced  in  humid  conditions  (see  Table  4). 
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Figure  30 .  Plots  of  Hydrazine  Concentration  Versus  Time  in  the  Teflon  Chamber 
(Four  Sequential  Runs  in  Humid  Air). 
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Figure  31.  Plots  of  Hydrazine  Concentration  Versus  Tine  in  the  Teflon  Chamber 
(Six  Sequential  Runs  in  Humid  Air). 
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TABLE  5.  HYDRAZINE  DECAY  IN  HUMID  AIR 


Date.  o£  Run  Half-Life  (hours) 

1.0 

1.5 

1.8 

3.3 

1.5 

1.5 

1.0 

1.3 
2.1 
3.1 
2.7 
3.0 

"Denotes  first  run  on  this  date,  etc. 

The  6-7  Jan  86  run  was  the  13th  to  be  made  in  the  5-mil  Teflon  chamber. 
The  previous  12  runs  had  little  conditioning  effect,  as  evidenced  by  the  1.5 
hour  half-life  of  the  6-7  Jan  run.  The  humidity  source  in  this  run  was  pure 
air  bubbled  through  distilled  water.  The  humidity  level  was  lower,  about  40 
percent,  and  the  temperature  of  the  room  was  unusually  low,  16  °C. 

The  four  sequential  experiments  conducted  on  21  and  22  Apr  86  were  made 
at  about  50  percent  relative  humidity.  There  was  an  increase  in  half-life 
from  the  first  run  made  on  21  Apr  to  the  second  run  made  on  that  same  day. 
However,  when  the  next  run  was  made  on  22  Apr,  the  half-life  decreased  by  50 
percent,  indicating  rapid  desorption.  The  final  run  on  22  Apr  showed  a 
slight  further  decrease  in  half-life.  As  in  the  high  humidity  run  conducted 
in  nitrogen,  there  was  an  approximately  fourfold  increase  in  ammonia 
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production  in  these  experiments,  compared  with  hydrazine  decay  in  dry  air. 
Figure  32  shows  fits  of  the  effusion-adsorption  version  to  the  first  two  runs 
in  this  sequence. 

The  final  six  sequential  experiments  conducted  on  10-12  Feb  87  were  run 
at  about  60  percent  humidity.  Humidity  variations  during  each  experiment 
were  less  than  2  percent,  except  for  the  final  two  runs,  where  the  humidity 
decreased  about  8  percent  over  the  course  of  the  run.  The  first  four 
experiments  showed  an  expected  increase  in  half-life  as  the  chamber  became 
conditioned  by  the  hydrazine. 

The  possibility  that  the  rapid  loss  of  hydrazine  in  both  humid  nitrogen 
and  humid  air  runs  could  be  due  to  association  of  hydrazine  with  water  in  the 
vapor  phase  to  form  hydrazine  hydrate  was  also  considered.  However,  infrared 
spectra  in  either  system  show  no  absorbance  features  which  cannot  be 
accounted  for  in  terms  of  hydrazine  or  water.  This  finding  is  supported  by 
an  early  study  of  Coulter  (Reference  23),  who  also  found  no  spectroscopic 
evidence  for  hydrazine  hydrate  formation.  Hence,  hydrazine  surface 
adsorption  is  most  likely  being  enhanced  under  humid  conditions.  This  being 
the  case,  the  effusion-adsorption  model  version  modified  for  conditioning 
closely  describes  what  is  actually  occurring. 

D.  CONCLUSIONS 

Experimental  and  modeling  results  show  that  gaseous  hydrazine  is  not 
oxidizing  at  ambient  conditions,  but  interacting  strongly  with  surfaces  in 
the  Teflon  chamber  environment.  The  model  described  in  this  report  fakes 
into  account  all  types  of  elementary  reactions  and  surface  interactions. 
Different  versions  of  the  model  have  been  used  to  fit  data  runs.  The  moat 
accurate  version  seems  to  be  the  effusion-adsorption  version  that  allows  for 
second-order  rather  than  first-order  adsorption.  The  sequential  data  sets 
show  a  conditioning  effect  aptly  described  by  second-order  adsorption. 
Binding  sites  become  almost  full  after  the  second  or  third  run  in  a 
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sequential  set  of  runs  in  dry  conditions.  Of  course,  this  process  is 
reversible,  so  that  the  chamber  can  be  deconditioned. 

Table  4  shows  results  of  the  effusion-adsorption  model  version  fits  to 
sequential  sets  in  dry  and  humid  conditions.  The  reproducibility  of  the 
fitting  parameters,  k^  and  [F,]o,  is  good  for  the  two  dry  runs  and  the 
two  humid  runs.  The  reproducibility  of  k  ^  is  not  as  good.  This  parameter 
may  vary  strongly  with  humidity  and  temperature.  This  seems  to  be  the  case, 
because  there  is  a  difference  in  the  values  of  k  ^  for  the  50  percent  and 
60  percent  humidity  runs.  The  parameter,  k  ,  varies  in  the  model  fits  for 

G 

humid  runs  as  well.  Again,  this  parameter  also  appears  to  be  sensitive  to 

environmental  conditions,  which  implies  that  k  does  not  just  represent 

6 

effusion,  but  steady-state  permeation,  which  no  doubt  is  somewhat  sensitive 
to  humidity. 

Interestingly,  k.,  k  and  k  increase  substantially  with 

X  A  ,  V 

hmidity.  That  k,  and  k  increase  implies  that  hydrazine  has  a  stronger 
affinity  for  a  wetted  surface.  However,  that  k_^  increase  shows  that 
adsorbed  hydrazine  comes  off  the  surface  sites  more  easily  in  humid 
conditions.  This  may  imply  that  with  the  wetted  surfaces  inside  the  chamber, 
hydrazine  ia  binding  more  to  the  water  molecules  rather  than  the  surface 
sites.  Whereas  the  hydrazine  has  a  strong  affinity  for  wetted  surfaces,  it 
apparently  binds  more  tightly  with  the  dry  surface.  This  can  be  described  by 
a  physical  process  in  which  hydrazine  is  becoming  trapped  in  the  structure  of 
the  surface,  requiring  more  activation  energy  for  adsorption.  In  terms  of 
energetics,  the  adsorption  and  desorption  activation  energies  are  lower  for 
wetted  surfaces,  as  illustrated  by  Figure  33.  This  makes  both  kj  and  k_j 
larger  for  humid  conditions.  That  k_1  ia  much  larger  in  humid  conditions 
shows  that  deconditioning  can  be  greatly  expedited  by  flushing  the  chamber 
with  humid  air.  This  is  reflected  in  Table  5  by  the  fact  that  allowing  the 
chamber  to  be  filled  with  humid  air  overnight  essentially  deconditiona  the 
chamber.  Otherwise,  using  dry  air  for  purging  would  require  a  few  weeks  of 
purging  to  decondition  the  chamber.  The  half-life  for  deconditioning  in  dry 
air  is  about  3  weeks,  whereas  in  60  percent  humid  air,  the  half-life  for 
deconditioning  is  2  hours.  As  might  be  expected,  [F1Jq  is  slightly 
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REACTION  COORDINATE 


Figure  33.  Diagram  Showing  the  Decrease  in  the  Activation  Energies  for 
Adsorption  and  Desorption  in  Wet  Conditions  (E  )  Compared  to 
Dry  Conditions  (E  ^).  A 
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less  in  humid  conditions,  implying,  of  course,  that  the  water  molecules  are 
occupying  the  sites  otherwise  available  for  the  hydrazine. 

The  fitted  values  of  [^]Q,  in  units  of  absorbance,  can  be  multiplied 

by  an  appropriate  conversion  factor  to  render  the  number  of  active  sites  per 

unit  area  of  Teflon  film,  provided  the  Teflon  film  is  the  most  important 

source  of  sites.  Given  the  following  values:  Teflon  surface  area  -  2.5 
2 

meter  ,  and  approximate  conversion  between  hydrazine  absorbance  and  number 
of  moles  -  3.8  x  10-5  mole'Au-*  (Au  =  absorbance  units)}  the  number  of 
sites  per  unit  area  can  be  calculated  to  be  .05  X  .  The  reciprocal  of 
this  is  20  X2.  For  humid  conditions,  this  number  increases  to  31  X2,  in 

o2  * 

good  agreement  with  the  value  of  23  A  found  by  Hayes  for  the  adsorption 
of  hydrazine  hydrate  on  certain  soil  components.  Of  course,  there  are  other 
surfaces  present,  such  as  the  mirrors,  which  also  seem  to  provide  adsorption 
sites. 

The  value  for  the  parameter,  k  [S] ,  from  the  permeation  model  version 

s 

does  increase  slightly  for  the  run  with  the  additional  2-mil  sheets  added 

(Table  4).  However,  k  [S]  is  significantly  larger  for  the  run  with  the 

s 

mirrors  added.  Although  the  permeation  version  is  not  an  accurate  model 

version,  the  fact  that  k  [S]  is  larger  is  significant.  Interactions  of 

s 

hydrazine  with  surfaces  other  than  Teflon  also  seem  likely  in  light  of  the 
fact  that  ammonia  is  produced.  It  is  unlikely  that  the  N-N  bond  rupture 
required  to  form  NH^  will  occur  either  in  the  gas  phase  or  on  Teflon 
surfaces.  Also,  the  fact  that  NH^  production  is  increased  in  humid  con¬ 
ditions  lends  support  that  catalysis  on  surfaces,  unlikely  to  be  Teflon,  is 
occurring,  since  initial  surface  adsorption  is  increased  in  humid 
conditions.  Moreover,  the  degree  of  NH3  production  is  insensitive  to 
oxygen  content,  giving  further  evidence  that  hydrazine  is  adsorbed,  or 
perhaps  even  chemisorbed,  onto  a  surface  which  is  able  to  rupture  the  N-N 
bond,  again  unlikely  to  be  Teflon. 
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Of  all  gases  studied,  hydrazine,  water,  and  0^  are  the  only  gases  which 
seem  to  interact  strongly  with  Teflon  film,  evidenced  by  rapid  changes  in  gas 
concentration.  Gases  chemically  similar  to  hydrazine,  such  as  methanol  and 
NH^,  surprisingly  do  not  interact  as  strongly.  Moreover,  the  molecular 
size  does  not  correlate  with  loss  rates  as  indicated  by  the  near  identical 
behavior  of  CH^  and  SF The  exact  chemical  or  physical  nature  of  the 
interactions  of  these  gases  with  Teflon  is  not  known  at  present.  However,  it 
is  apparent  that  surface  interactions  and  physical  loss  processes  are  the 
most  important  loss  processes  of  gaseous  hydrazine  in  Teflon  chambers  at 
ambient  conditions. 
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APPENDIX  A 


LISTING  OF  MACRO  PROGRAM  (AND  SUB-MACRO  PROGRAMS)  TO  AUTOMATICALLY 
ACQUIRE,  PROCESS,  STORE,  PRINT  AND  PLOT  FT-IR  DATA 


This  program  is  called  TS3.  It  was  written  in  the  Nicolet  MACRO  language 
by  using  the  TED  editor  and  then  compiling  with  the  Nicolet  utility  program 
MACCRT.  This  procedure  allows  remarks  to  be  included  in  the  program  for 
documentation  and  user  assistance. 

The  main  MACRO  program  (TS3)  calls  several  SUB-MACRO  programs.  Most  of 
these  were  written  and  compiled  as  subsections  of  TS3.  Each  program  ends  with 
the  statement  END.  Each  new  program  begins  with  followed  by  the  program 
name. 


\MACRO  TO  ACQUIRE,  PRINT,  AND  PLOT  SPECTRA 
YAUTQMAT I CALLY 


ITS3 

BFN 

DFN 

NSD 

OMD 

OMD 

SELECT  BEGINNING  AND  ENDING  FREQUENCIES  FOR  PPK,  REGION  1 

VI 0 

VII 

THRaO.005 

OMD 

SELECT  BEGINNING  AND  ENDING  FREQUENCIES  FOR  PPK,  REGION  2 
VI 2 
VI 3 
OMD 

- 4  INITIAL,  SHORTER  INTERVAL  SPECTRA  WILL  BE  ACQUIRED 

OMD 

ENTER  THE  DELAY  TIME  (MIN)  BETWEEN  THESE  INITIAL  SPECTRA 

VF1 

OMD 

- THE  REMAINING  SPECTRA  WILL  B°  ACQUIRED  AT  A  LONGER  INTERVAL 

TSX 

END 

\\\\\\\ 


JTSX 

OMD 

ENTER  THE  DELAY 

VFO 

OMD 


\MACRO  TSX;  (TS3  WAS  TOO  LONG,  TSX  IS 
\THE  2ND  PART) 

TIME  (MIN)  BETWEEN  REMAINING  SPECTRA 
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OMD 

PRESS  ••RETURN"  TO  BEGIN 
OMD 


PAU 

CSN=0  \SET  THE  DUMMY  VARIABLE  CSN  TO  ZERO 

ISS  \RUN  THE  SUB-MACRO  ISS 

CSN=0  \SET  THE  DUMMY  VARIABLE  CSN  TO  ZERO 

TSS  \RUN  THE  SUB-MACRO  TSS 

END 


IISS  \SUB-MACRO  TO  ACQUIRE,  PROCESS  AND  PLOT  THE 

\FIRST  4  SPECTRA 
\ 

CSN=17  \SET  DUMMY  VARIABLE  FOR  UPCOMING  CMR  COMMAND 

FOR  BBB=1  TIL  4 

AQ2  \SUB-MACRO  TO  ACQUIRE,  PROCESS  AND  PRINT  DATA 

CMP=CSN  \SET  CONDITIONAL  MACRO  PARAMETER  FOR  EACH  PLOT 

CMR 

DFN=DFN+1 

CSN-GSN+1 

DLI  \SUB-MACRO  TO  EXECUTE  SOFTWARE  TIME  DELAY 

NXT  BBB 

END 


•DLI 

OMD 


\SUB-MACRO  TO  GIVE  TIME  DELAY  FOR  INITIAL  4 
\ SPECTRA  IN  "TS3" 


OMD 


OMD 

OMD 


DELAY  IN  PROGRESS 


OMD 

DCL=204445  \SET  NEW  DISPLAY  COLORS 

DCX  \ CHANGE  DISPLAY  TO  NEW  COLORS 

FOR  111=1  TIL  VF1 

FOR  AAA=0  TIL  60000  \LOOP  TO  DELAY  FOR  1  MIN 

NXT  AAA 
NXT  III 

DCL=210414  \RETURN  COLORS  TO  DEFAULT  VALUES 

DCX  \EX£CUTE  COLOR  CHANGE 

END 


1M17  \SUB-MACRO  TO  PLOT  UPPER  LEFT  SPECTRUM 

\PART  OF  TS3 
\ 


XSP=VF2 

XEP=VF3 

PL0=100 

ZPN 

XPN=1 

YPN=5 . 5 

PEN 

XSL=3 

YSL=4 

PAG=NO 

TIT=NO 

PLS 

XPN=.5 

YPN=4.5 

PEN 

PON 

PRN  DFN 
STD 
POF 
END 


!M18 


XSP=VF2 

X£P=VF3 

XPN=>3.5 

YPN=-4.5 

PEN 

PLS 

XPN=.5 
YPNa4 . 5 
PEN 
PON 

PRN  DFN 
STD 
POF 
ZPN 

XPN=-4,5 

PEN 

ZPN 

END 


NSUB-MACRO  TO  PLOT  UPPER  RIGHT  SPECTRUM 
\PART  OF  TS3 


IM19  NSUB-MACRO  TO  PLOT  LOWER  LEFT  SPECTRIN 

\PART  OF  TS3 
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XSP-VF2 

XEP=VT3 

PLS 

XPN= . 5 
YPN=4 . 5 
PEN 
PON 

PRN  DFN 

STD 

POF 

XPN-3 . 5 
YPN=~4 . 5 
PEN 
END 


JM20 


XSP=VF2 

XEP=VF3 

PLS 

XPN=.5 
YPN=4 . 5 
PEN 
PON 

PRN  DFN 

STD 

POF 

ZPN 

XPN»3 

PEN 

2PN 

END 


\SUD-MACR0  TO  PLOT  LOWER  RIGHT  SPECTRUM 
\PART  OF  TS3 


1TSS  XSUB-MACRO  TO  COMPLETE  DATA  ACQUISITION 

\AN0  HARD  COPT  OUTPUT  FOR  MTS3“ 

\ 

FOR  AAA=1  TIL  100  \LOOP  TO  ACQUIRE  AND  PROCESS  100  FILES 
CSN»17 

FOR  CCC»1  TIL  4  \LOOP  TO  ACQUIRE,  PROCESS,  PRINT  AND 

\PLOT  4  SPECTRA 

AQ2  XSUB-MACRO  TO  ACQUIRE  AND  PRINT  DATA 

CMPaCSN 

CMR 

DFNaDFN+1 

CSN=GSN+1 

DL2  XSUB-MACRO  TO  EXECUTE  THE  DESIRED  DELAY 

\BETWEEN  SETS  OF  SCANS 

NXI  CCC 
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NXT  AAA 
END 


!DL2  \SUB-MACRO  TO  EXECUTE  THE  DESIRED  DELAY 

NBETWEEN  SETS  OF  SCANS 
\PART  OF  TS2 
OMD 

OMD 


OMD 

OMD 


DELAY  IN  PROGRESS 


OMD 


DCL=204445 

DCX 

FOR  T1I=1  TIL  VFO 

FOR  AAA=0  TIL  60000 

NXT  AAA 

NXT  III 

DCL=.210414 

DCX 

END 

J! 


The  SUB-MACRO  AQ2  was  written  and  compiled  separately  and  is  listed  below. 


•AQ2  \SUB-MACR0  TO  ACQUIRE,  DO  A  PPK  AND  PRINT  PPK  VALUES 

\ 

TKM=>4  \SET  TEM  TO  TOGGLE  PRINTER  ON 

NPR  \TURN  ON  PRINTER 

SCD 

PRN  DFN 

STD 

SFN=2 

RFN=DFN 

OFNaDFN 

MOS 

RAS 

ABS 

DFN=SFN 

XSP=VIO 

XEPaVIl 

PPK 

XSP=VI2 
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XEP=VI3 

PPK 

OHD 


DFN=RFN 

XSP=4000 

XEP=400 

YSP=-.l 

YEP=1 . 1 

DSS 

TEM=4 

RPR 

END 

1 1 
•  « 


\SET  TEM  TO  TOGGLE  PRIKTER  OFF 
\TURR  PRINTER  OFF 
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APPENDIX  B 


LISTING  OF  MACRO  PROGRAM  (AND  SUB-MACRO  PROGRAMS)  TO  CALCULATE 
BASELINE-CORRECTED  ABSORBANCE  VALUES  FOR  HYDRAZINE,  METHANE 
AND  SULFUR  HEXAFLUORIDE 

This  program  is  called  MPK.  It  was  written  in  the  Nicolet  MACRO  language 
by  using  the  TED  editor  and  then  compiling  with  the  Nicolet  utility  program 
MACCRT.  This  procedure  allows  remarks  to  be  included  in  the  program  for 
documentation  and  user  assistance. 

The  main  MACRO  program  (MPK)  calls  several  conditional  MACRO  programs 
which  were  written  and  compiled  as  subsections  of  MPK.  Each  program  ends  with 
the  statement  END.  Each  new  program  begins  with  followed  by  the  program 
name.  The  final  program  in  the  series  is  terminated  with  "l  l"  after  its  END 
statement. 

This  listing  also  includes  the  FORTRAN  program  BASE1,  which  is  called  by 
the  MACRO  during  the  baseline  correction  routine. 


:MPK  VMACRO  TO  MEASURE  NEW  ABSORBANCE  BY  CHOOSING  A 

\ BASELINE  FOR  EACH  ANALYTICAL  BAND  FOR  EACH  SPECIES 
\ 


OMD 

CHOOSE  SPECIES  TO  BE  ANALYZED: 
OMD 

OMD 


HYDRAZINE  a 
METHANE  « 


SF6 


CMP 

OMD 

ENTER  BACKGROUND  FILE  NO.  FOR  THIS  ANALYSIS 

BFN 

OMD 

ENTER  BEGINNING  FILE  TO  BE  ANALYZED 

DFN 

OMD 

ENTER  NO.  OF  FILES  TO  BE  ANALYZED 
QIT 

QIT=QIT+DFH 

SRToDFN 


30 

31 

32 


SFN=2 

FOR  XXX- SET  TIL  QIT 

RAD 

ABD 

ASD 

OMD 

FILE  CURRENTLY  BEING  ANALYZED  IS: 
PRN  DFN 


DSD 

CMR 

DFN=DFN+1 
NXT  XXX 
E1ID 


IM30  \CONDITIONAL  MACRO,  PART  OF  MPK 

\THIS  MACRO  DETERMINES  THE  BASELINE  AND  NET 
\ ABSORBANCE  FOR  HYDRAZINE  AT  -958  CM-1 
OMD 

BASE1  FRSQ=DBB=  1170  CM-1  (HYDRAZINE) 

OMD 

BASE2  FSEQ=RBB=  860  CM-1  (  M  ) 

OMD 

BASE  I  ABSORBANCEsVFl ,  BASE?.  ABSORBANCE-TBB 

XSP=1180 

XEP=1160 

DSD 

TEM-6  \TURN  ON  PRINTER 

NPR 

PAU 

TEMsS  \REC0RD  ABSORBANCE  AND  FREQ.  VALUES  WITH  CTRL-W 

NPR 

PAU 

TEM«7  \TURN  OFF  CURSOR 

NPR 

PAU 

DBB=PXA  \STORE  VALUE  OF  PXA  IN  DUMMY  VARIABLE  DBB 

VF1=PYA  \STORE  VALUE  OF  HA  IN  DUMMY  VARIABLE  VF1 

TEM=A 
NPR 

XSP=860 

XEP=840 

DSD 

TEM=>6  \TURN  ON  CURSOR 

NPR 

PAU 

TEM=8  \ RECORD  ABSORBANCE  AND  CM-1 

NPR 

PAU 

TEM=7  \TURN  OFF  CURSOR 

NPR 

PAU 

RBB=PXA  \STORE  NEW  PXA  VALUE  IN  DUMMY  VARIABLE 

TBB=PYA 

XSP=958  \SBT  X-AXIS  VALUES  TO  BRACKET 

XEPo957  \  THE  HYDRAZINE  ANALYTICAL  PEAK 

THR=0 
PPK 


RBB 
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S&£s*XA 

RAN=PYA 


GFN 

BASHl.FTH 

SIB 

ra? 

TEM=4 

NPR 

END 


\STORE  NEW  PXA  VALUE  IN  DUMMY  VARIABLE  PXA 
\STORE  NEW  PYA  VALUE  IN  DUMMY  VARIABLE  PYA 

\CALL  TILJ  FORTRAN  BASELINE  CORRECTION  PROGRAM 

\RUN  THE  FORTRAN  PROGRAM 


!M31  \CONDITIONAL  MACRO,  PART  OF  MPK 

\THIS  MACRO  DETERMINES  THE  BASELINE  AND  NET 
UBSORBANCE  FOR  METHANE  AT  -1306  CM-1 

OMD 

BASE1  FREQ=DBB=  1320  CM-1  (METHANE) 

OMD 

BASE 2  FSEQ=RBB=  1310  CM-1  (  "  ) 

OMD 

BASE1  ABSORBANCEaVFl ,  BASE2  ABSORBANCE«TBB 

XSPsl330 

XEP=1310 

DSD 

TSN»«  \TUSN  ON  POTIgg. 

NPR 

PAU 

TEM=8  \WG0S3>  ABSORBANCE  AND  FREQ.  VALUES  WITH  CTRL-W 

NPR 

PAU 

TE8*7  \t\m  on  CURSOR 

NPR 

PAU 

08WXA  \STORE  VALUE  OF  PXA  IN  DUMHY  VARIABLE  D8B 

WlsPYA  \ STORE  VALUE  OF  PYA  IN  DUMHY  VARIABLE  VY1 

NPR 

XSP=1295 

XEPel275 

DSD 

TBH=6  \TUSN  ON  CURSOR 

NPR 

PAU 

TEK=8  \ RECORD  ABSORBANCE  AND  CH-1 

NPR 

PAU 

TEMa7  \TURN  OFF  CURSOR 

NPR 

PAU 

EBB=PXA  \STORE  NEW  PXA  VALUE  IN  DUMMY  VARIABLE  EBB 

TBBsPYA  \STORE  NEW  PYA  VALUE  IN  DUMMY  VARIABLE  IBB 
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XSP-1306  \SET  X-AXIS  VALUES  TO  BRACKET 

XEP=1305  \  THE  METHANE  ANALYTICAL  PEAK 

THR=0 
PPK 

OGN=PXA  \STORE  NEW  PXA  VALUE  IN  DUMMY  VARIABLE  PXA 

RAN=PYA  \STOSE  NEW  PYA  VALUE  IN  DUMMY  VARIABLE  PYA 

GFN 

BASE1.FTN  \CALL  THE  FORTRAN  BASELINE  CORRECTION  PROGRAM 

STD 

FRN  \RUN  THE  FORTRAN  PROGRAM 

TEM-4 

NPR 

END 


:M32  \CONDITIONAL  MACRO,  PART  OF  MPK 

\THIS  MACRO  REPLACED  THE  BASELINE/NET  ABSORBANCE 
\METH0D  (AS  IN  M31  AND  M32)  WITH  THE  STANDARD 
\NIC0LET  INTEGRATION  ROUTINE,  SET  TO  COVER  THE 
\948  CM- I  BAND  OF  SF6. 

OMD 

INTEGRATION  MACRO 
SFN=2 

ENTER  BACKGROUND  AND  DESTINATION  FILES 

BFN 

DFN 

OFN^DFN 

HOS 

RAS 

ABS 

FXFalOOO 

LXF»900 

BAS=YS 

SMS 

TEM»4 

NPS 

PRN  DFN 
STD 

PRN  FCS 
TEMs»4 
NPR 
END 


The  final  Hating  ia  the  FORTRAN  program  BASB1.FTN. 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

C  PROGRAM  BASEl!  TAKES  ABSORBANCE  AND  FREQUENCY  VALUES 


cioo  ooo  ooonoooo  o  o  o  o  o  n  o  o  o 


FROM  FT-IR  AND  DOES  A  LINEAR  INTERPOLATION 
TO  DETERMINE  NET  ABSORBANCE  FOR  AN 
ANALYTICAL  PEAK. 


REAL  DBB,VF1,OGN,RAN,RBB,TBB,VFO,VF2,VF3,VF4,FOV 

RETRIEVE  BASELINE  AND  ABSORBANCE  VALUES  FROM  FT-IR 
USING  THE  RVAL  FUNCTION 

DBB=RVAL( 6216,1) 

VF1=RVAL( 6136,1) 

OGN=RVAL< 6249,1) 

2AN=RVAL( 6247,1) 

RBB=RVAL( 6218,1) 

TBB=RVAL (6214,1) 

CALCULATE  NET  ABSORBANCE  VALUE 

NET  ABSORBANCE  (FOV)  =  PEAK  ABSORBANCE  (RAN)  -  EXTRAPOLATED 

BASELINE  ABSORBANCE  (VF4) 

VF4  =  (VF1-TBB)  X  ( OGN-RBB) / ( DBB-RBB ) 

VFOsABS(TBB-VFl) 

VF3aABS ( DBB-RBB ) 

VF2»ABS ( DB3-0GN) 

VF4-VF2/VF3 

VP4WF4+VF1 

F0V*EAN-VF4 


PRINT  OUT  RESULTS  IN  FT-IR 
WRITE  (2,40)  FOV 

40  FORMAT  (40(,-,),/»’  V  FEAK  ABSo* #F9.3,/,40( '-'),/) 

RETURN  TO  FT- IE 

CALL  IE2XIT 
END 
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APPENDIX  C 


LISTING  OF  MACRO  PROGRAM  TO  CALCULATE  THE  BASELINE-CORRECTED  ABSORBANCE 
VALUES  FOR  HYDRAZINE  (957.3  CM-1  BAND),  INCLUDING  THE  FORTRAN  PROGRAM 

GET2.FTN 

This  program  is  called  BC1.  It  was  written  in  the  Nicolet  MACRO  language 
by  using  the  TED  editor  and  then  compiling  with  the  Nicolet  utility  program 
MACCRT.  This  procedure  allows  remarks  to  be  included  in  the  program  for 
documentation  and  user  assistance. 

The  main  MACRO  program  (BC1)  calls  the  FORTRAN  program  GET2.FTN,  which,  in 
turn  calls  a  FORTRAN  subroutine  called  FGETOF.  These  programs  are  listed,  in 
sequence,  below. 


IBCl  \MACR0  TO  COMPUTE  BASELINE-CORRECTED  ABSORBANCE 

WALUES  FOR  HYDRAZINE  (957.3  CM-1  BAND) 

\ - USED  WITH  SINGLE  BEAM  SPECTRA 

BFN 

OMD 

ENTER  BEGINNING  AND  ENDING  DFN'S 

SRT 

QIT 

DFN=SRT 

SFN=2 

FOR  AAAcSRT  TIL  QIT 

TEM=4 

NPR 

PRN  DFN 
STD 
TEM=4 
NPR 

OFN=DFN 

MOS 

HAS 

ABS 

DFN=SFN 

XSP»958 

XEP*957 

TllR=-.5 

PPK 

XSP=4000 

XEP=0 

TEM=t4 

NPR 

GFN 

GET2 . FTN 

FRN 

TEM=4 
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NPR 

DFN=OFN 
DFN=DFN+1 
NXT  AAA 
END 


The  FORTRAN  program  GET2.FTN  is  listed  below. 

CCCCCCCCGCCCCCCCGGCCGCCCCCCCGCCCCGCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCGCCCCC 

G 

C 

C  PROGRAM  GET2.ASC 

C 

C  THIS  PROGRAM  RETRIEVES  ABSORBANCE  DATA  FOR  HYDRAZINE 

C  FROM  THE  DESTINATION  FILE  AND  STORES  IT  IN  A  REAL 

C  FORTRAN  ARRAY,  RDATA(I),  FOR  USE  IN  BASELINE 
C  CORRECTIONS.  IT  IS  CALLED  FROM  THE  FT-IR  MACRO  BC1. 

C 

C  THE  BASELINE  CORRECTION  ROUTINE  IS  AN  EXTENSION  OF  THE 

C  LINEAR  CORRECTION  GIVEN  IN  MFT-I8  SPECTRAL  LINES",  VOL. 

C  II,  NO.  1,  PP.  8-9. 

C 

C  (D.A.S.  22  DEC  86) 

C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C234567 

PROGRAM  GET2 

INTEGER  IHEAD(512),ISTP,IENP,NP 
REAL  RDATA<  23000) ,BP1 , LBP , PXA , PYA , CABS ,XDIF,XDIF2 
1,YDIF,FAC1,FAC2,BP,SP 

COMMON  /  BUFBL  /  IDA(23000),IBUF(512) 

RPMAXb23000 

C 

C  GET  THE  DESIRED  REGION  OF  THE  DFN  INTO  A  REAL  FORTRAN  ARRAY 
C 

CALL  FGETOF( I HEAD , I STP , IENP ,NP,NPHAX , FXAX, RDATA) 

C 

C  RETRIEVE  THE  FILE  HEADER  INFORMATION  AND  CHECK  FOR 
C  CORRECT  SPECTRUM  FORMAT  AND  ARRAY  SIZE. 

IF  (IHEAD(IO).NE.O)  THEN 
PRINT  100 

100  FORMAT  (•  V SPECTRUM  NOT  ABSORBANCE',/) 

CALL  IREXIT 

ELSEIF  (NP.GT.NPMAX)  THEN 
PRINT  101,  NPMAX 
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uoo 


101  FORMAT  (’  ','N0.  OF  PTS.  REQUESTED  IS  MORE  THAN  \I7,/) 
CALL  IREXIT 
ENDIF 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COMPUTE  THE  BASELINE  CORRECTED  ABSORBANCE  VALUE  FOR  THE 
ANALYTICAL  PEAK  OF  THE  DESIRED  SUBSTANCE: 

THE  STARTING  AND  ENDING  POINTS  FOR  DETERMINING  THE 
BASELINE  ARE  FOUND  BY  AVERAGING  THE  ABSORBANCE  VALUES 
IN  A  SPECIFIC  RANGE  FOR  THE  BEGINNING  AND  ENDING 
POINTS  OF  THE  BASELINE.  THE  LOCATION  OF  THESE 
POINTS  IN  THE  SCRATCH  FILE  IS  DETERMINED  FROM  THE 
FACT  THAT  THE  DATA  FILE  IS  44  SECTORS  LONG  WITH  512 
POINTS  (OR  WORDS)  PER  SECTOR,  I.E.,  22528  POINTS 
(THE  VALUE  OF  FSZ  FOR  A  1.0  CM-1  SPECTRUM).  HOWEVER, 
ONLY  16384  POINTS  ARE  NECESSARY  TO  DEFINE  A  1.0  CM-1 
SPECTRUM.  IN  ADDITION,  WITH  A  SAMPLE  SPACING  (SSP) 

OF  2,  THE  FREQUENCY  RANGE  COVERED  BY  THE  DATA  IS 
ZERO  TO  1/2  THE  LASER  FREQUENCY  (0-7899.1  CM-1). 

THUS,  WITH  16384  DATA  POINTS  TO  COVER  7899.1  CM-1, 
THERE  ARE  (7899.1/16384)  0.482  CM-1  PER  DATA  POINT 
IN  THE  FT-IR  SCRATCH  FILE.  NOW,  TO  DETERMINE  THE 
DATA  POINT  VALUES  FOR  A  GIVEN  WAVENUMBER  RANGE,  THE 
ABOVE  RELATIONSHIP  IS  USED.  FOR  EXAMPLE,  IN  THE 
HYDRAZINE  BASELINE,  THE  FIRST  BASELINE  POINT  WILL 
BE  THE  AVERAGE  VALUE  OF  THE  POINTS  FROM  1160  TO  1180 
CM-1,  OR  2406  TO  2488  IN  THB  DATA  FILE.  SO  THB 
AVERAGE  OF  THE  POINTS  (43)  ION  THIS  PART  OF  THE  DATA 
FILE  WILL  BE  DETERMINED  BY  ADDING  THEM  UP  AND 
DIVIDING  BY  43. 


AVERAGE  THE  POINTS  FOR  THE  FIRST  BASELINE  VALUE 


SP=0 
BP  1=0 
COUNT=0 

XFAC=7899. 1/16384 
DO  200  1=2406,2448 
WAVE=XFAC*(I+1) 
SP=SP+RDATA(I) 
COUNToCOUNT+1 

200  CONTINUE 

BP1=SP/C0UNT 


AVERAGE  THE  POINTS  FOR  THE  SECOND  BASELINE  VALUE 


EP=0 

LBPaO 


94 


o  o  o  o  o  n  on  non  onoooo  ooooo 


COUNT=0 

DO  201  1=1742,1784 
WAVE1=XFAC*  ( 1+1) 

EP=EP+RDATA(I) 

C0UNT=C0UNT+1 

201  CONTINUE 

LBP=EP/ COUNT 

RETRIEVE  THE  FT-IR  VALUES  FOR  FREQUENCY  AND  ABSORBANCE 
FROM  TH  PEAK  PICK  ROUTING.  (THEY  ARE  STORED  AS  FLOATING 
POINT  VALUES  IN  PXA  AND  PYA) 

PXA=RVAL( 6112,1) 

PYA=RVAL(6116,1) 

CALCULATE  BASELINE  CORRECTION.  THERE  ARE  TWO  POSSIBLE 
CASES;  (A)  POSITIVE  SLOPING  BASELINE  OR  (B)  NEGATIVE 
SLOPING  BASELINE. 


CASE  A  -  POSITIVE  SLOPING  BASELINE 
(I.E.,  BF1  <  LBP) 

SKIP  TO  CASE  (B)  IF  BASELINE  SLOPE  IS  NEGATIVE 
IF(BPZ.GT.LBP)  GOTO  3000 
OTHERWISE,  PROCEED  WITH  THE  CORRECTION 


YDIFaLBP-BPl 
XDIFs.1170-850 
KDIF2=1170-PXA 
FAC1*XDIF2/XDIF 
FAC2=FAC1*YDIF 
FAC2=FAC2+BP1 
GOTO  3500 

CASE  B  -  NEGATIVE  SLOPING  BASELINE 
<1.8. ,  BP1  >  LBP) 


3000  CONTINUE 

YDIF=BP1-LBP 

XDIF=1170-850 

XDIF2=PXA-1170 

FAC1=*XDIF2/XDIF 

FAC2«FAC1*YDIF 

FAC2=FAC2tLBP 

c 


3500  CONTINUE 
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c 

C  APPLY  CORRECTION  TO  ABSORBANCE  VALUE  FROM  PPK 
C 

CABS=PYA-FAC2 

C 

C  PRINT  OUT  CORRECTED  ABSORBANCE  VALUE 

C 

WRITE  (2,550)  CABS 

550  FORMAT(40( *-'),/, ‘  ' , ‘PEAK  ABS=  • ,F9.3,/,40( •-* ),/) 
C 

C  RETURN  TO  FTIR 

C 

CALL  IREXIT 
END 


The  final  listing  is  the  FORTRAN  subroutine  FGETOF, 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c  SUBROUTINE  FGETOF 
C 

C  Gets  a  file  status  block  IKBAD  and  piece  of  spectrum  EDATA 
C  from  FT-IR  data,  using  XSP,  XEP,  OFN,  and  FS2S  (in  sectors). 

C 

C  DWV,  Nicolet  -  Dec.  1981  (revised  Jan,  1982,  DWV) 

C 


C 

C 

C 

C 

C 

C 


Revision  attached  to  NICOLET  letter  dated 
28  Aug  86  (completed  2  Dec  86),  D.A.S. 

Revised  to  have  a  form  compatible  with 
the  original  FGSTOF  program  (Nicolet 
FORTRAN  77,  pp.  61-62),  D.A.S.  9  Dec  86 


C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C234567 

SUBROUTINE  FGETOF  (IKRAD,X$TP,IENP,NP,NFHAX,  FXAX,  AND  RDAXA) 

COMMON  /  BUFBL  /  IDA( 23000) 

REAL  RDATA( 23000) 

INTEGER  IHRAD(512),DFN,PSZS 
FLASBRal5798.2 
C 

C  Use  IRVAL  to  load  necessary  values  from  FTIR  parameters 
C 

IXSP  *  IRVAL(6165,1) 

IXBP  «=  IRVAL(6166,1) 

DFN  =  IRVAL( 6162,1) 
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cioo  ana  onn  n  o  o 


FSZS  =  IRVAL( 6145,1) 

INODE=  IRVAL( 13003,0) 

Sow  get  file  status  block  of  DFN  file  and  call  it  IHKAD 

CALL  IRTISK( I HEAD , 352 , 87+(DFN+l)*FSZS , INODE) 

NPTS  =  IHEAD(15)*128 
IFGF  =  IHEAD(8) 

FBXP  =  2** ( 19-IHEAD ( 6) ) 

FXAX  =  FLASER/ (NPTS*IHEAD( 17 ) ) 


Now  calculate  ISTP  and  IENP  for  interferogram  or  wavenumber  case 

IF  (IFGF.EQ.O)  THEN 
FXAX  =  1. 

IF  (IXSP.LE.XXEP)  THEN 
ISTP  =  ISXP 
IENP  «  IXEP 
ELSE 

ISTP  a  IXEP 
IENP  a  IXSP 
ENDIF 

ELSBIF  (ISEP.GT.IXEP)  THBN 
ISTP  a  IXEP/FXAX+G , 5 
IENP  a  IXSP/FXAX+0.5 
ELSE 

ISTP  «  IXSP/FXAX+0.5 
IENP  «  IXEP/FXAX+O . 5 
ENDIF 

HP  IENP-ISTP+1 

Check  for  a  reasonable  number  of  points  in  HP 

IF  (NP.GT.NPMAX)  THEN 
PRINT  901,NPKAX,NP 
HP  =  HPKAX 

ELSEIF  (KP.LT.5)  THEN 
PRIHT  902, HP 
BHDIP 

How  read  data  from  scratch  DFN  into  array  IDA 
IN  a  1 

HBONE  a  ISTP/512 
NSKIP  =  ISTP-NDONE*512 

IF  (HSKIP.GT.O)  THEN 

CALL  IRTISK(IDA, 512, 88+DFN+FSZ+NDONE, INODE) 

DO  10  XaKSKIP+1,512 
IDA(IN)  »  XDA(I) 
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o  o  o  o  o  o  o  o  o 


IN  =  IN+1 


10  CONTINUE 

NDONE  =  NDONE+l 
ENDIF 

NPREK  =  NP-IN+1 
IF  (NPREM.GT.O)  THEN 

CALL  IRTISK(IDA(IN) , NPREM , 88+DFN*FSZ+ND0NE , itfODE ) 

ENDIF 

la  FTIR,  (-524288)  is  the  "official  ridiculous  value"  which 
is  used  in  blanking  spectral  areas.  Since  FORTRAN  doesn't 
and  shouldn't  recognize  this  value,  a  trick  is  used  to  create 
this  value  in  IRID  in  order  to  excise  it  from  RMTA. 

I RID  =  -524287-1 

Nov  float  IDA  into  RDATA,  checking  for  IRID  and  zeroing. 

DO  50  1=1, NP 
IDAS  =  IDA(I) 

IF  ( IDAS. EQ, IRID)  IDAS  =  0 
RDATA(I)  a  IDAS/FEXP 
50  CONTINUE 

RETURN 

901  FORMAT  (IX, 'Only  dimensioned  to  ',16,'  points,  not  ,16/) 

902  FORMAT  (IX, 'Why  are  there  only' ,16, 'points  requested?'/) 

END 
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APPENDIX  D 


THE  PERMEATION  MODEL  VERSION  AND  THE  BASIC  SUBROUTINE  USED  IN  THE 
MARQUARDT  NONLINEAR  LEAST  SQUARES  ROUTINE 

The  permeation  model  allows  for  permeation  through  the  chamber  walls  as 
the  only  kinetic  process.  Hence,  from  Equation  (11),  the  model  reduces  to: 

d[G]/dt  =  -ka[S][G]  +  k^IG-SJ 

d[G-S1]/dt  =  kg[S][G]  -  k_fl[G->Sj]  +  k^SlUG-S^  -  [G-S^) 

For  i  =  2  to  N-l  (if  N  A  3); 

dlG-Sjl/dt  =  kptSKlG-S^}  +  EG-S1+13  -  216-Sj]) 

d(G-SR)/dt  »  -  (G-Sh]>  -  k_a[G-Sjj) . 

Since  it  ia  assumed  that  (S)  >>  (G) ,  (S)  is  essentially  constant  and  the  only 
three  parameters  are  lcg CSJ ,  k  g,  and  k^lSj .  N  ia  a  predetermined 
integer.  Although  in  theory  N  is  quite  large,  in  practice  it  is  normally  set 
to  2  for  expediency  in  calcuiati-jns.  Making  N  larger  does  alter  the 
parameters,  but  does  not  severely  affect  the  goodness  of  the  fit.  For  N  »  2, 
the  dependent  variables  are  [G] ,  (6-S^] ,  and  tG-S^l . 

The  actual  differential  equations  which  have  been  used  in  the  routine  are: 

dL^/dt  =  cj(g  -  Li)  -  cj(Li  -  L2) 

For  i  ®  2  to  N{ 

dLj/dt  a  C3(Lj_i  +  L*+i  -  2Lj) 


dLu+i/dt  =  C2<lN  -  lN+1>* 

in  which  cx  =  kg[S],  c2  "  k-s*  c3  3  VS1,  and  the  1,8  are  given 
in  the  following  expressions: 


CK) 


EG]  =  g  -  Li 


For  1  *  l  to  Hj 

[G-Si]  =  Li  -  Li+i« 

These  equations  are  generated  from  the  following  analysis: 


G  +  Si  .4.  —  ,  G-Si 

g  -  Li  k_3  Li 


For  i  =  2  to  N; 


G—S i _ 1  +  Si 
Li_i  -  Li 


kD 

+  -  —  G-$i  +  si-l 

kD  M 


k_s 

G-Sj|  . . ■">  G  (outside  chamber)  +  Sg. 

lN  -  LN+1  Lh+I 

Note  that  EG],  which  is  proportional  to  the  absorbance  according  to  Beer's 
Law,  la  the  variable  being  monitored.  This  analysis  allows  for  a  direct 
monitoring  of  how  much  of  G  is  lost  from  the  chamber  system.  It  can  easily 
be  seen  that: 

R 

g  -  Lu+i  *  CGI  ♦  £  EG-Sj]. 

i=l 

L^j  is  the  concentration,  in  moles  per  unit  chamber  volume  of  gas  lost 
from  the  system.  It  cannot  be  interpreted  as  the  concentration  outside  the 
system,  since  all  concentration  terms  must  have  r  common  reference  point,  In 
this  case  the  chamber  volume. 

The  following  is  the  subroutine,  which  starts  on  1100  of  the  program, 
used  to  perform  the  numeric  integration: 

1100  If  I  >  1,  then  1125 


loo 


1X01  Dt  =  .01 

1102  Nm  =  2 

1103  Dim  L<Km+l) 

1104  For  E  *  1  to  Km  +  1 

1105  L(E)  =  0 

1106  Next  E 

1107  T  a  0 

1108  T  =  T  +  Dt 

1109  L(l)  *  LCD  +  C(1)*(G  -  L(l))*Dt  -  C(2)*(L<1)  -  L(2))*Dt 

1110  For  E  »  2  to  Nm 

1111  1(B)  a  L(E)  +  C(3)*(L(E-1)  +  i»(E+l)  -  2*L(E))*Dt 

1112  Next  E 

1113  LCNm+1)  *  t<Nm+l)  +  C(2)*(L(Nm)  -  L(No4-i))*Dt 

1114  U  T  <  0.1  Then  UOS 

1115  If  ABS(G  -  LCD  -  "initial  absorbance  value")  <  0.001  then  1118 

1116  G  s  L(l)  +  "initial  absorbance  value" 

1117  Go  to  1104 

1118  II  a  G  -  L(l) 

1119  T  «  0 

1120  Go  to  1140 

1125  T  «  T  ♦  Dt 

1126  1.(1)  *  I»(l)  4  C(1)*(G  -  L(l))*Dt  -  C(2)*(L(1)  -  L<2)*Dt 

1127  For  E  »  2  to  Km 

1128  L(E)  •  L(E)  +  C(3)*a<B-.l)  +  t(E4l)  -  2*L(E))*Dt 

1129  Next  E 

1130  L(Ka4l)  =  L(Ka+l)  +  C(2)*(L(Hm)  -  L(N»+l))*Dt 

1131  If  T  <  X<I,1)  then  1125 

1132  Y1  s  G  -  L(l) 

1140  Return 

The  terms  are  defined  as: 


Km  -  number  of  pore  layers 

Y1  -  dependent  variable,  in  this  case  the  absorbance  readings 
T  -  time 

Dt  -  differential  time  element,  At 


101 


X(I,1)  -  independent  variable,  in  this  case  the  Ith  time  reading.  The 
"1"  refers  to  this  being  the  first,  and  in  this  case  the  only,  independent 
variable. 

The  subroutine  operates  in  the  following  manner.  Initially  I  is  set  to 
one,  so  the  program  first  goes  through  lines  1101  to  1120.  These  lines 
define  the  parameters  Dt,  Nm,  and  T,  and  set  first  conditions.  In  the 
experimental  runs,  the  gas  is  mixed  in  the  chamber  between  5  and  10  minutes 
before  the  first  absorbance  reading  is  taken.  This  is  the  mixing  time. 
Hence,  enough  time  has  elapsed  to  allow  for  some  surface  adsorption,  and 
consequently  G  in  lines  1109  and  1126  is  not  the  initial  absorbance  reading, 
but  is  slightly  larger.  Lines  1108  through  1117  allow  integration  to  occur 
for  the  mixing  time,  which  is  usually  set  at  0.1  hour  for  convenience.  Since 
G  is  not  known  exactly  at  the  beginning  of  the  program,  it  is  generally  set 
to  the  initial  absorbance  reading  or  slightly  larger  before  program 
execution.  Line  1116  readjusts  G  and  the  numeric  integration  for  the  mixing 
time  is  repeated.  This  is  done  sufficient  times  so  as  to  make  the  quantity, 
G  -  L(l) ,  within  a  certain  tolerance  (.001  absorbance  unit)  of  the  initial 
absorbance  value.  Once  this  condition  has  been  met,  T  is  reset  to  zero,  I  is 
set  to  2  and  the  program  goes  to  line  1125  of  the  subroutine.  Here  the 
integration  begins  with  T  =  0  and  goes  to  the  full  time  of  the  run.  The  main 
part  of  the  program  then  calculates  the  sum  of  squares,  readjusts  the 
parameters,  and  goes  through  the  entire  integration  subroutine  again.  The 
intent  of  the  program  is  of  course  to  minimize  the  sum  of  squares. 
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APPEND TX  E 


ANALYSIS  OF  SEVERAL  KINETIC  SYSTEMS 


kl 

1,  System:  G - ►Products 

k2 

2G - ►  Products 


Differential  equation:  d[G]/dt  =  -(k2[G]  +  k]_)[G], 


Integration: 

I  d[G]/{[G](k2[G]  +  k!»  =  -  /  dt, 


[G]  «  lkigexp(-k1t)]/{k2g[l  -  expC-k^)]  +  k^ 


kl 

System;  G  + 

A  — 

B 

g-x 

a-x 

X 

k2 

G  + 

B 

- ►  Product 

g-x~y 

x--y 

y 

Differential  equations: 

d[A]/dt  =  -kitGHA],  or  dx/dt  =  k^(g  -  x  -  y)(a  -  x) 
dIB]/dt  =  -k2[G)[Bj  +  kilGJU},  or 

dx/dt  -  dy/dt  o  -k2(g  -  x  -  y)(x  -  y)  +  ki(g  -  x  -  y)(a  -  x) 

Therefore,  dy/dt  =  k2(g  -  x  -  y)<x  -  y).  Y  can  be  solved  in  terms  of  x 
as  follows: 

dy/dt/dx/dt  =  dy/dx  a  k^(x  -  yj/tk^a  -  x)J. 


10.1 


Let  a  =  k2^ki*  If  dy/dx  is  written  as  dy/dx  =  -ay/(a  -  x)  +  ax/(a  - 
x)  an  integrating  factor  can  be  easily  introduced  to  aid  in  the 
integration.  This  factor  is; 

exp[a/dx/(a  -  x)3  =  (a  -  x)-a. 

Hence,  (a  -  x)  a  dy/dx  +  (a  -  x)-a-1ay  =  (a  -  x)-a-1ax,  or  d[y(a  -  x)-a]  = 
(a  -  x)-a-1axdx.  The  expression,  (a  -  x)-a_1  axdx,  can  be  integrated 
by  parts  to  yield: 

aJ*(a-x)_Q-1xdx  =  x(a  -  x)~a  -  a  -  x)~adx. 

The  overall  expression  can  now  be  integrated  for  the  following  two  cases: 
Case  1:  a  »  1; 

/d[y/(a  -  x)]  »  x/( a  -  x)  -  /dx/(a  -  x)  +  G, 

in  which  C  is  the  integration  constant.  At  t  »  0,  x  »  y  =»  o.  Hence,  C 
=  -lua  and  y  is: 

y  *>  x  +  (a  ~  x)ln(l  -  x/a). 

[G]  and  [B]  can  now  be  put  in  terms  of  [A] : 

[G]  =  g  -  2(a  -  [A])  -  [A3ln(tA)/a), 

[BJ  *  tA]ln(a/[A] ) . 

Case  2:  a  ji  1; 

Jd[y(a  -  x)’°)  »  x/(a  -  x)°  -  /(a  -  x)~0  dx  +  C. 

In  this  case,  C  =  a1’0/^  -  1).  Therefore  y  is: 
y  a  [a  -  ax  -  a1_°(a  -  x)a]/(l  -  a). 
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[G]  and  [B]  can  be  shown  to  be: 

[5]  =  g  -  a[l  -  <[A]/a)a  +  (1  -  2a) (1  -  [A]/a)]/(l  -  a) 

[B]  =  [([A]/a)a-l  -  1] [A]/(l  -  a). 

The  differential  equations  in  time  can  be  analytically  integrated  easily 
when  a  =  1/2.  In  this  case,  the  equation  for  d[A]/dt  is: 

d[A]/dt  =  -k1[AJ{g  -  2a[l  -  ([Al/a)1^)}. 

If  g  =»  2a,  d[A]/dt  =  -2k!a[A]([A]/a)1/f2  and  integration  yields: 

i'lAl  At 

/  a~1/2  [A]“3/2  d[A]  =  -2k!  /  dt, 

[A]  *  a/[(l  +  Iqat)2] 

[6]  =  2a(l  -  k!at)/(l  +  kjat) 

[Bj  =  2k!a2t/[(l  +  kiat)2]. 

1/2 

If  g  *  2a,  making  the  substitution,  8  »  ([A] /a)  ,  yields  the 

following  simplified  equation: 

de  =  -k!8[g  -  2a(l  -  0)]dt/2. 

Integration  within  limits  yields: 

rfl  /t 

g  -  2&<1  -  8))>  .  -(kx/2)  j  dt, 

8  =>  (2a  -  g)/[2a  -  gf(t)], 

in  which  f(t)  =>  expi.(g'2  -  a)k!t].  The  concentrations  car  be  shown  to  be: 

U3  «  a{(2a  -  g)/[2a  -  gf(t)]>2 
IG)  »  g  +  2ag[f<t)  -  l]/[2a  -  gf(t)] 

IB]  *  2ag(2a  -  g)[l  -  f(t)]/[2a  -  gf(t)]2. 
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3.  System:  G  +  A  Zj: 


kl 


g-x  a-x  k_i  x 


Differential  equation: 


d[G]/dt  =  -k].[G][A]  +  k_i [B] ,  or  dx/dt  =  kx(g  -  x)(a  -  x)  -  k_jx. 


The  integral  within  limits  is: 


-  x)(a  -  x)  -  k_ix] 


.  /"at. 


This  integral  can  be  found  in  mathematics  handbooks.  The  integration 
yields: 

i* 

(l/v'q)ln[(2x  -  a  -  g  -  k_x^l  ->/5)/(2x  -  a  -  g  -  k_x/kx  «  t, 

|Q 

in  which  q  =  (a  +  g  +  k_x/kx)*  -  4ag,  which  is  always  positive  unless  a 
■=  g  and  k_x  =  0*  In  this  case  it  is  zero.  X  can  be  shown  to  be: 


x  =  Z1Z2[exp<k1>/qt)  -  l]/{2[Zx  -  Z2exp(klN/qt)]>, 

in  which  Zxa^-a-g  -  k_x/kx  *nd  Z2  =  -(\/q  +  a  +  g  +  k>x/kx) . 
The  concentrations  are: 


(B]  *  2ag[f(t)  -  l)/[Zx  -  Z2f(t)] 

[A]  -  a{l  -  2g[f(t)  -  13/[2x  -  Z2f(t))> 
[GJ  =  g{l  -  2a(f(t)  -  l]/[Zx  -  Z2f<t)]}, 

in  which  f(t)  =>  exp(kxv^t). 

kl  . 

4.  System:  G  +  A  B  +  C 

g-x  a-x  k_x  x  x 
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Differential  equation: 


d[G]/dt  =  -k]_[G][A3  +  k_].[BHC],  or  dx/dt  =  k^(g  -  x)(a  -  x)  -  k_]X2. 

Integration  within  limits  yields: 

^dx/[(ki  -  k_i)x2  -  (a  +  g)k]X  +  k^ag]  =  t. 

If  =  k_i,  x  can  be  shown  to  be: 

x  =  ag{l  -  exp[-(a  +  g)kit]>/(a  +  g). 

The  concentrations  are: 

[B]  =  [C]  =  ag[l  -  f(t)]/(a  +  g) 

[6]  =  g[g  +  af(t)]/(a  +  g> 

[A]  =  a[a  4  gf(t)]/(a  +  g), 

in  which  f(t)  =  exp[-(a  4  g)k^tl.  If  k^  /  k  the  form  of  the 
integral  ia  identical  to  that  of  system  3.  Therefore,  the  concentrations 
can  be  shown  to  be: 

[BJ  a  [C]  =  2ag[f(t)  -  Z2f(t)] 

[G]  «  g{l  -  2a(f(t)  -  l]/[Zi  -  Z2f(t)]> 

[k]  »  a{l  -  2g{f(t)  -  l)/[Zx  -  Z2f(t)J), 

2  r~ 

in  which  f< t)  =  exp(klNftt),  q  «  (a  -  g)  +  Aagk^/k^  Z^  »  -  a  -  g, 

and  Z2  o  -  (V?  4  a  +  g). 

kl  . 

5.  System:  G  +  A  ZjZ ZZZfc  B  Assumption:  a  >>  g 

g-x  a  k_i  x 

k2 

A  4  B  3=fc  C 
a  x-y  k_2  y 
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Differential  equations: 


d[G]/dt  =  -kia[G]  +  or  dx/dt  =  k^a(g  -  x)  -  k_i(x  -  y) 

d[C]/dt  =  k2a[B]  -  k_2[C],  or  dy/dt  =  k2a(x  -  y)  -  k„3y 

X  and  y  can  be  put  in  terms  of  one  another  as  follows: 

dx/dt/dy/dt  =  dx/dy  = 

[kia(g  -  x)  -  k_i(x  -  y)]/[k2a(x  -  y)  -  k_2yl * 

The  differential  equation  can  be  solved  by  making  the  following 
substitutions: 

k2ax  -  (k2a  +  k_2)y  ^(x  +  c2)  +  c3(y  +  C4) 

kiag  -  (kxa  +  k^Jx  +  K.&  =  c5(x  +  c2)  +  c6(y  +  C4), 

in  which  the  constants  can  be  shown  to  be: 

ci  o  k2a 

c2  =  -k1ag(k2a  +  k.2)/[ki«(k2a  +  k»2)  ♦  *-1^-21 
c3  =  -k2a  -  k_2 

04  o  -kik2a2g/[kjia(k2a  +  K.2)  +  k-ik-g] 

CS  a  -kia  -  k_i 
c6  «  k_i. 

The  differential  equation  reduces  to  : 

dy’/y’  =  (K2a  -  K2  -  l)ds/{-K2s2  +  [K2  -  1  -  o(l  +  K^la  +  a}, 

in  which  y*  =  y  +  c^,  0=  k_j/k_2S  =  k^a/k_^,  c 

k2a/k  2,  x'  s  ay’,  and  x'  =  x  +  c2«  The  variable,  a,  has  been 

introduced  to  allow  for  variable  separation  in  the  equation.  It  can  be 
shown  that  the  numerator  of  the  right-hand  side  of  the  equation  above  is 
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never  the  derivative  of  the  denominator.  The  integral  of  this  expression 
can  be  found  in  mathematics  handbooks,  and  is: 


lny*  =  -(l/2)ln(-K2s2  +  es  +  a)  +  [(9/2  -  K2  -  l)/yq]ln[(2K2s  -  e 
+  ^/q)/(2K2s  -  0  ~yjq) ]  +  C, 

In  which  e=K2-l-a(l  +  K^) ,  q  =  a2  +  4K2a  ,  and  C  is  the 
integration  constant.  Att  =  0,x  =  y  =  0.  Solving  for  C  and 
substituting  back  into  the  expression  yields: 

y  =  K^gtl  -  f<s)3 /[K1(K2  +  1)  +  1], 

in  which  f(s)  «  <[K2(1  +  1/K2)2  -  9(1  +  1/K2>  -  o]/(K2b2  -  es  - 

a  )}1/2{(2K2s  -  0  +>/q)(2  +  2^-9  -<y(q)/[(2  +  2K2  -  0  + 

•^q)(2K2s  -  o-  and  4>  =  (9/2  -  K2  -  l)/^q.  Noting  that  x  = 

s(y  +  cA)  -  c2>  the  concentrations  can  be  put  in  terms  of  s  as: 

[B]  »  KiK28(1  -  a)f(8)/[K1(K2  +  1)  +  1]  +  K1g/[Kji(K2  +  1)  +  1] 

[C]  a  KxKjjgll  -  f(s)l/[K1(K2  +  1)  +  1] 

[GJ  «  K1K2gsf(s)/[K1(K2  +  1)  +  1)  +  g/tKi(K2  +  1)  +  1]. 

Relative  values  for  these  concentrations  can  be  calculated  by  inputting 
the  appropriate  values  for  s.  The  upper  value  of  s  is  the  value  at  t  »  0. 
At  t  a  0,  s  a  1  +  1/K2>  At  t  o  infinity,  the  concentrations  are  all  at 
their  equilibrium  values.  The  equilibrium  concentrations  can  be  shown  to 
be; 

[B3e  «  Klg/[Ki(R2  +  1)  +  1] 

[C]e  =  K^g/lKid^  +  1)  +  1] 

[G]e  «  g/[Ki(K2  ♦  1)  +  1J, 

in  which  the  subscript  "eM  denotes  equilibrium.  Equilibrium  conditions 
are  met  where  f(s)  a  0.  Since  $  is  always  negative,  f(s)  <=  0  where 
2K2s-  0  •¥yjh  =  0,  Hence,  sfl  =»  (9  \/5)/2K2 ,  The  range  for  s  is: 


1  +  1/K2  2  S  2  (e  +Vq)/(2K2). 

2 

It  is  true  that  the  term,  K5s  -  es  -  a  ,  i3  2ero  where  s  =  s  . 

4  6 

Even  though  this  term  appears  in  the  denominator,  the  term,  2K2s  -  e 
is  in  the  numerator,  and  f(s)  does  approach  zero  as  s  approaches  sg. 

The  differential  equations  in  time  can  be  solved  by  assuming  a  gen¬ 
eral  exponential  form  of  the  solution  to  the  differential  equations  in 
time.  This  is  done  in  several  kinetics  textbooks  and  will  not  be  done 
here. 
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APPENDIX  F 


ADAPTATION  OF  THE  MODEL  FOR  DESCRIBING  EFFUSION  AND  PERMEATION 

INTO  THE  CHAMBER 


For  a  gas  coming  into  the  chamber,  one  can  write  the  following  series  of 
steps: 


ke 

G(outside)  G(inside) 

g  ke  Le 


G(outside)  +  Sjj  T  .T  G-Sr 
8  k-s  lN+1 

For  i  =  1  to  N  -  1; 

kp 

G-Sn+1-i  ♦  sN-i  + - £  G'SN-i  +  SN+1.i 

LN+2-i  -  kD  LN+l-i 

LN+l-i 


k_g 

G-Sx  +=£  G<  inside)  +  Si 
l2  ~  ^1  ka  L1 


This  model  assumes  an  infinite  outside  reservoir  of  the  gas.  The  concentra¬ 
tions  are: 


[G]  63  Le  +  Lx 


For  i  »  1  the  Nj 

[G-S^]  a  Li+x  —  Lx. 


The  differential  equations  are: 


dLe/dt  a  kg(g  -  Le  -  Lx) 
ffltqf+l/dt  “  k8tS]g  -  k_a(LR+x  -  Ljj) 
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For  i  =  1  to  N  -  1; 

dLN+l-i/dt  "  kDtsK%+2-i  +  LN-i  '  2LN+l-i> 
dLj/dt  =  k_s(L2  -  ^i)  -  ka(Le  +  i*i)» 
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APPENDIX  G 


ADAPTATION  OF  THE  MODEL  FOR  DESCRIBING  PERMEATION  THROUGH 
CHAMBER  WALLS  AND  THROUGH  PERMEABLE  SURFACES  WITHIN  THE  CHAMBER 

For  a  gas  permeating  through  the  surface  of  a  material  placed  inside  the 
chamber,  the  processes  are  described  as: 

_  ^s  _ 

G(inside)  +  j-  ?  G-Si 

g  -  Li  K.a  Li 


G-Si+i  +  ^ 
^i+1 


k_B 

inside)  +  Sjj, 
ka  g  -  Li  + 

%+l 

in  which  the  overhead  bar  refers  to  the  permeable  surface  inside  the 
chamber.  For  a  gas  permeating  through  the  chamber  vails,  the  processes  are: 

ka 

G(inaide)  +  Sj^'- G-Sj 
B  -  M  k.8  Li 

£fi+l  “  L1 

For  i  *  l  to  H  -  lj 

kp  ^ 

G-Si  ♦  S1+1^Zr*:G-Si+i  ♦  Si 
Lt  -  k0  L1+1 

Li+1 

k-s 

G-Sfl  ’  G(outslde)  +  Sjj . 

lN  ~  lN+1 

lN+1 


%*■! 


For  i  a  1  to  N  -  1; 

*D 

G-Si  +  SU1 
Li  -  kD 

Ci+1 
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The  concentrations  are: 


[G]  =  g  -  Li  +  LJj+i  -  L]_ 

For  i  =  1  to  N; 

[G-Sj]  =  Li  -  Li+i 

For  i  9  l  to  N; 

[G-SiJ  =  Li  -  Li+i» 

The  differential  equations  describing  the  gas  permeation  through  the  inside 
permeable  surfaces  are: 

d[G]/dt  =>  -dLx/dt  +  dLg+]/dt  -  dL1/dt 

=  -(2kg[S]  +  ka[S])(g  -  L1  +  L^+1  -  Lj) 

+  Wh  k*  i>  ♦  k-.<Li  - 

d[G-B1)/dt  ■=  dLj/dt  -  dE  /dt 

=  -  h  ♦  £g+1  -  y  -  -  y  -  ^[SkEj  -  ♦  i3> 


For  i  a  2  to  K  -  1; 

dlG-Sjl/dt  a  dEj/dt  -  dL1+1/dt  *  kjjts] [Lt-1  -  Hlt  -  tul)  -  li<r2} 
d[G-S-]/dt  a  dLg/dt  *  dL^+]/dt 

=  V5»<»  -  h  ♦  -  y  -  *-.<%  -  %♦!> 

-  ysHLjj^  -  2Eg  +  %*!>• 


The  Independent  variables  can  be  combined  algebraically  in  pairs  and  the 
number  of  differential  equations  reduced  to  yield: 

dL1/dt  «  2k8[5](g  -  Lx  ~  L^>  -  k  ^  -  t2) 
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For  i  =  2  to  q  +  6  -  1  (if  N  2  5); 
dL^/dt  =  kD[S3(Li_1  +  L*1+1  -  2^) 

-  V°> [E,+i,-i  - 2(1  +  ‘>Ei+.i  <lf  s  2  3)- 

in  which,  in  general,  Lj  =  Lj  -  L^+2i,  1  *  N/2,  and  6  =  1/2  if  N  is 
odd  and  zero  if  N  is  even.  The  equations  for  the  permeation  through  the 
chamber  valla  are: 

dLj/dt  =  ka(S](g  -  L1  -  Lr)  -  k_s(L1  -  L2) 

For  i  =  2  to  H; 

dl^/dt  »  kD[S](Li_1  +  L1+1  -  2l»i) 

dLN+l/dt  °  k-a(LN  ’  Hl+l5, 

Since  L’^  =  [G]  0  *  "  h  ’  ^1’ 


APPENDIX  H 


FIT  OF  HYDRAZINE  TT-UAY  DATA  FROM  STUDIES  AT  THE  UNIVERSITY  OF  LEEDS 

Hydrazine  decay  studies  in  Pyrex  containers  were  conducted  at  the 
University  o;T  Leeds  (Reference  10).  No  attempt  was  made  to  model  the  kinetic 
data  by  the  researcher  at  the  university.  Some  of  the  data  from  the 
University  of  Leeds  studies  have  ueen  modeled  in  this  work.  Two  models  have 
been  used  to  test  the  concepts  outlined  in  this  report  against  the  data:  the 
adsorption-vapor  oxidation  model  and  the  surface  catalysis  model.  The 
adsorption-vapor  oxidation  model  is: 


N2H4  +  S  N2H4-S 

k_s 


N2H4  +  02 


"►  i'2  +  2H20, 


in  which  the  rate  constants  have  their  usual  meaning.  The  surface  catalysis 
model  ic; 


N2H4  +  S  «  N2H4-S 

*-a 


n2r4-s  +  o2 


N2  +  2H20  +  S 


If  the  variable,  x,  represents  the  extent  of  the  first  reaction  in  terms  of 
the  partial  pressure,  and  y  represents  the  extent  of  the  second  reaction,  the 
differential  equations  'or  the  adsorption-vapor  oxidation  model  can  be 
written  as: 


dx/dt  =  kslS](Pfl2H4  -  x  -  y)  -  k_sx 


dy/dt  =  kp(Pj52H4  -  x  -  y)<Po2  ~  y) 


116 


in  which  the  Fs  represent  initial  partial  pressures.  The  data  was  taken  as 
the  total  pressure  vs  time.  In  terms  of  initial  partial  pressures  and  x  and 
y,  the  total  pressure  is  given  from  this  model  as: 

pt  =  pn2h4  +  po2  -  x  +  y 

The  differential  equations  for  the  surface  catalysis  model  are: 
dx/dt  »  k8[SKpN2H4  -  x)  -  k_s(x  -  y) 
dy/dt  =  ke(x  -  y)(Po2  -  y) 

The  total  pressure  for  this  model  is  given  by: 


PT  =  pN2H4  +  Pq2  “  x  +  2y 


If  the  approximation  is  made  that  P 


°2 


y  **  PA  end  k  **  0, 
02  -s 


then  the  surface 


catalysis  model  can  be  analytically  Integrated  to  yield: 


y  =  PN2H4-S  +  pi:2H4  -  kpiiBifxvi-l'WQtb  -  k‘) 

+  lk,pH2H4'<kp  -  k’>  ~  pM2H4-Sl«*P<-kpt) 

in  which  k*  *>  katS},  k£  =  kpp02»  411(1  PN2H4-S  14  partial  pressure  equivalent 
of  the  initial  surface  concentration  of  N2H4-S. 

Both  models  were  fitted  to  representative  data.  Figure  H-l  shows  a  data 
set  with  model  fits.  As  can  be  seen,  the  surface  catalysis  model  fits  the 
data  much  better.  Results  of  the  fit  tor  the  surface  catalysis  yield:  k'  a 

.0‘4( .002)  minute-1,  kppQ2  a  .0325(.008)  minute-1,  and  Pn2H4-S  “  .18(.05)  mm  Hg. 
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1‘igure  H-i.  Model  Fits  to  University  of  Leeds  Hydrazine  Decay  Data. 

Conditions:  PM  H  a  10  am  Hg,  Pn  =  140  ram  Hg,  Temperature 

n2m4  u?  -1 

26  °C,  Surface-co-Volume  Ratio  (Pyrex  Vessel)  »  1.38  cm  . 


